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INTRODUCTION 


K. Schwarzschild’s classical work on absorption and dififiision in 
the sun’s atmosphere, and the continuation of these investiga- 
tions through A. S. Eddington, J. H. Jeans, E. A. Milne and others, 
have rendered the theory of radiative equilibrium a definite 
chapter of mathematical astrophysics. The problems connected 
with this theory are of the following type. (Jaseous material 
filling a certain part of space is subject to given incident radia- 
tion.* What is the distribution of light and of temperature in the 
medium, when in radiative equilibrium? In the outer layers of 
a star, where the curvature can be neglected, i.e. where the 
material can be considered as stratified in parallel planes, those 
problems appear in their simplest mathematical form. The 
astrophysicists mostly contented themselves with an a})proxi- 
mate solution of these ])roblems. Owing to a certain inherent 
beauty, however, they aroused also the interest of the ligorous 
mathematician. It is the purpose of this tract to attempt a 
coherent representation of all that has been achieved in the 
direction of a rigorous treatment of those standard problems. 

The rigorous solution of a i)hysical problem naturally pre- 
supposes an exact formulation of the physical assumptions im- 
plied in the problem. Has then the rigorous solution more than 
merely mathematical significance, when these assumptions j)rove 
right only with a limited degree of approximation? I’he answer 
to this question is not always in the negative. The {Schuster- 
Schwarzschild model of a purely scattering atmosphere, or the 
Schwarzschild model of a purely absorbing gray atmosphere, 
both in radiative equilibrium, represent for instance typical 
standard models which play a similar role in the theory of the 
stellar atmospheres as the ‘intermediate orbits’ in celestial 
mechanics. The Milne model, being of special importance as a 

* More generally, on the boundary the radiation satisfies given conditions. 
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simple limit case (infinite optical depth) of Schwarzschild’s model, 
may well be compared with Hill’s periodic orbit in the lunar 
theory. It is well known what the detailed study of an inter- 
mediate orbit means for the computation of the actual one. 
A more detailed and rigorous treatment of the above standard 
models of radiative equilibrium should thus deserve more than 
mere mathematical interest. 

Besides the known results (theory of Milne’s standard model, 
infinitely deep slab with no radiation incident on its surface) the 
reader will also find more or less detailed discussions of other 
models. Milne’s model with radiation incident on the boundary 
(insolation of a planetary atmosjjherc, reflexion effect in close 
binary stars), the Schuster-Schwarzschild model with an 
arbitrary law of scattering, narrower limits for the solution of 
the Schuster-Schwarzschild standard problem, Milne’s model of 
a planetary nebula, and other problems. 

From the mathematical point of view, the main feature of the 
above models is that they are governed by integral equations, 
which are linear in the standard cases (pure absorption and gray 
material, j)ure scattering). Its important property, namely the 
positivity of the kernel, has already been recognized by Schwarz- 
schild as fundamental for the rigorous discussion. It may be 
mentioned that the reading of the book re([uires no special know- 
ledge of the theory of integral equations, since the positivity of 
the kernel allows all problems to be treated in an elementary 
way. Only Chapter iv presupposes some knowledge of the 
Fourier integral. 

It is a pleasure to the author to thank Prof. Milne for many 
helpful remarks and the Syndics of the University Press for 
accepting this book as a Tract. 

E. H. 

.Voi’. 1, 1933 



CHAPTER I 


FUNDAMENTAL PRINCIPLES AND PROBLEMS 

§1. TRANSMISSION OF RADIATION THROUGH 
ABSORBING AND SCATTERING MATERIAL 

The outer layers of a star are constantly exposed to an enormous 
flow of radiant energy coming from the deep interior. In working 
itself through the outer layers, a part of the radiation is absorbed 
and scattered. The scattered radiation will immediately be re- 
distributed in the different directions, the frequency being 
unchanged. On the other hand, a particle heated up by the 
absorption of radiation re-emits temperature radiation of all 
frequencies. 

As a first approximation, the state of matter and radiation 
may be considered stationary. Therefore we must first find 
the conditions for a steady state. The intensity of radiation 
changes in a definite way along the ray. A part will be lost through 
absorption and scattering. On the other hand, the intensity gains 
again since the particles on the ray emit scattered and tempera- 
ture radiation in the direction of the ray. This law of transfer holds 
separately for each frequency. 

Combined with the condition for a steady state the energy prin- 
ciple yields another equation. The net loss of radiant energy 
within an element of volume (radiation of all frequencies) equals 
the net gain of heat energy through conduction and convection 
plus the radiant energy liberated from other (sub-atomic) sources. 
Energy of the latter kind must play a role in the interior of the 
stars in order to maintain their radiation over cosmical ranges of 
time. 

When energy is transferred by radiation only, the net flow of 
radiation through each volume element vanishes. This state of 
affairs is usually called radiative equilibrium (in the strict sense). 
Modern astrophysics claims the outer layers of a star to be 
approximately in radiative equilibrium. In the stellar interior 
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the energy liberated is no more negligible though being small in 
comparison to the radiant energy itself. In spite of this fact it 
has become customary to speak of radiative equilibrium (in the 
wider sense) in this case. 

For the determination of the temperature distribution, one 
must know the laws relating the temperature of matter to the 
radiation. When matter and radiation are in an enclosure, 
Kirchhoff’s and Plank’s laws hold. In the stars, however, the 
temperature varies from layer to layer, so that those laws cannot 
be strictly valid. In the interior, however, the radiation is prac- 
tically enclosed, and even in the outer layers those laws have 
proved to be a very good approximation. 

In the steady state all quantities describing the state are in- 
dependent of time. The radiation field is described by the intensity 
of radiation j 

as a function of the point P, the direction r and the frequency j^.* 
Its physical meaning is the following. Let da be an element of 
surface about P with the (sensed) normal n, and let dco be an 
infinitesimal bundle of directions containing r (surface element 


on the direction sphere). If we write 

dq~ I cos (it, r) \da (1) 

for the projection of da upon a plane perpendicular to r, the 
expression I^(P,r)dqdajdv (2) 


represents the radiant energy of the spectral interval dv which, 
per unit time, flows through da and spreads out in the solid angle 
doj. Here is meant, the bundle of all rays through all points of da 
and in all directions contained in dco, dq is evidently the perpen- 
dicular cross-section at P of the bundle. It is convenient to have 


a special name for the quantity 

I^(P,r)dqdv. (3) 

We call it the strength at P of the parallel bundle of the spectral 
interval dv, of direction r and with the cross-section dq. 

We set dadv^^^^^ = dadvlI^(P,r)co^(n,r)dw, (4) 

* The degree of polarization is not analysed here. 
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where the integral is extended through all directions r, and where 
n is a given direction. Let da be an orientated surface element 
about P and with n as its positive normal. The quantity 
represents then the net flux of radiant energy (of the spectral 
interval dv) through da and per unit time. Energy flowing through 
da from its negative towards its positive side (cos > 0) is here 
reckoned positive; negative, however, in the opposite directions 
r(cos < 0). Let x, y, z be three orthogonal directions, and let us 
write a, jS, y for the three direction cosines of the normal n with 
respect to x, y, z. The equation 

cos {n, r) = a cos (r, x) + jS cos (r, y) + y cos (r, z) 

shows that we have 

« = cos (w, x) + cos (/i, y) + cos (yi, z) , 

i.e. that is the component in the direction n of a vector 5^. 
This vector briefly called the net jlux at P of the 

v-radiation. 

In a vector held we have the Gauss integral identity 

j = (6) 

where V is a volume, S its surface, and n the outward normal of 
the surface element da. 

For the total radiation of all frequencies we set 

l = r I, dv, 5v= f (0) 

Jo Jo 

We now come to the quantities that describe the interaction 
between matter and radiation. Let rj^(P) be tlie ynass coefficient 
of emission at P and within dv. The quantity 

T]^{P)dmdojdv (7) 

represents then the radiation emitted, per unit time, by the mass 
element dm at P, within dv and within the solid angle doj. In 
analogy to the above we find it convenient to call 


7]^{P)d7ndv 


( 8 ) 
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the strength of a parallel bundle emitted by dm within dv. For 
the total emission of all frequencies we set 

yi(P) = ^\AP)dv. (9) 

Let, furthermore, (P) denote the mass coefficient of absorption 
for the v-radiation at P. Along a short path ds, the amount 
poLyl^ds of the intensity is absorbed, p being the density of 
matter at P. Outside of astrophysics it is more customary to 
introduce the linear coefficient of absorption, pa^ . 

The mass coefficient of scattering (P) is defined in an analogous 
way. Scattering weakens the intensity by the amount pa^l^ds 
along ds. pa^ is the linear coefficient of scattering. The radiation 
lost by scattering at P will immediately be redistributed among 
the different direction r' issuing from P. This distribution is 
described by a law of scattering 

j yAP-,r,r'), (10) 


where r denotes the direction of the incident ray. We have, of 

course, r /n nj / < /n\ 

( 11 ) 

The law of scattering is supposed to have the reciprocity property 

yp (P; r, r') = (P; r', r) (11') 

and ^ '*') = yA^> ^ (11") 


where —r denotes the direction opposite to r, i.e. the same 
amount of the scattered radiation is sent into r and into — r. 


These conditions are fulfilled under very general assumptions 
about the material. The simplest case y^= 1 (uniform scattering) 
has hitherto found the chief attention of the astrophysicists. It is, 
however, necessary to consider more general laws too, for instance 


Rayleigh’s law 


yv = y (^ ^') = -4 {1 + ^os2 (r, /)}. 


It is seen to fulfil (11') and ( 1 1"), as does also every law for which 
is a function of | cos (r, r') | . 

All quantities (except the net flux) are non-negative, 
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§2. ABSORPTION AND SCATTERINO OF THE MASS 
ELKMENT. THE EQUATION OF TRANSFER 


We now compute the v-radiation absorbed and scattered within 
a mass element dm at a point P. For this purpose let us first con- 
sider the parallel bundle of all rays through dm which have a 
given direction r. It is convenient to imagine dm, subdivided into 
many thin columns of direction r. Let ds be the length, dq the 
cross-section of such a column. The strength (P, r) dqdv of the 
parallel bundle of direction r through that column, i.e. through 
dq, is, in consequence of absorption within the column, weakened 
by the amount pcL^I^dqdvds. Since pdqds is the mass of the 
column, we infer by summation over all columns that 

I^dmdv 

is the amount that the strength of the bundle considered above 
loses on account of absorption within dm. By strength of that 
bundle is meant here the sum of the strengths of all the above 
partial bundles, once taken when entering dm, the other time 
taken when leaving dm. oi^I^,dmdv is the difference of these two 
sums. Now, energy is obtained by directional integration of the 
strength. Integration over all directions r of the above loss of 
strength shows thus that the whole di/-radiation travelling 


through dm per unit time loses the amount 

dmdvaL^{P)lI^(P,r)do} ( 12 ) 

in consequence of absorption. 

In the same way it is seen that the amount 

dmdvGy,{P)lIy(P,r)du) (13) 


is scattered within dm. We now show how this scattered radiation 

is redistributed among the different directions r\ According to 

the above considerations, the radiation flowing through dm per 

unit time and spreading out within a solid angle dco (containing r) 

loses the amount ^ j j t fn \ 

E = dmdva^(P) I^(P, r) 


through scattering within dm. This scattered amount E will be 
distributed over the other directions r' according to the law of 
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scattering. The part falling within a cone doj' about the direction 
r' is therefore ^ 


Summation over all cones tZoi, i.e. integration through all incident 
directions r, gives, therefore, 

dmdoj'dv'^''^P J I, (P, r) y, (P; r, r') dco (U) 

as the part of the scattered radiation (13) which is redistributed 
within the cone dco'. If we now integrate this over all directions 
r\ we obtain, according to (11), again (13) as it should be. 

We note for later f)urposes that the j)arallel bundle of scattered 
radiation emitted by dm into the direction r has the strength (the 
notation for incident and scattered ray is interchanged) 

dmdv J (P, r') y, ( P. r', r) doj’ (14') 

Equation of tra7isfer. Let us consider a light ray of direction r 
through P, and let P' be a nearby ])oint lying in the direction r 
from P, ds = PP'. We construct a thin cylinder about ds with the 
two bases at P and P', and with the cross-section dq. The strength 

I,{P',r)dqdv, (15) 

of the parallel bundle of direction r through dq consists then of 
three different parts. Firstly, the strength at P of the same bundle 
weakened by absor})tion and scattering in the cylinder, 

{1 - p (a„ + a^) ds} (P, r) dqdv; 

secondly, the strength (8) of the parallel bundle emitted by the 
cylinder in direction r, with 

d7n = pdqds; 


thirdly, the strength (14') of the bundle of scattered radiation 
sent by the cylinder into the direction r. On dividing the equation 
obtained in this way through pdqdsdv, we get Schwarzschikrs 
fundamental equation of transfer 


1 

p{^v + Ov) 


dlAP^r) 

ds 


= J,{P,r)-I^P,r). 


(16) 
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together with 

(«, + aj J, (P, r) = 7,, (P) + '"•’IP J/, ( P, /) y, ( P; r', r) dco'. 

(!«') 

The derivative on the left-hand side of (1(5) is the directional 
derivative with respect to P, taken in the direction r of the ray. 
On introducing three orthogonal axes x, y, z, P=(x, y, 2 :), and on 

setting a = cos(?%x’), ^ = cos(r,?/), y = cos(r, 2 ) (l^l^O 

for the direction cosines of r, we obtain 

dl dl dl 

^ 

The cpiantity J^(P,r) introduced by (16') may, according to 
Schwarzschild, be called the Ergiebigkeit, It is an average of the 
two quantities . ^ 


having thus the dimension of an intensity. The Ergiebigkeit 
depends, in general, upon the direction r. In the case of uniform 
scattering (y^=l) or of isotropic radiation (the intensity is in- 
dependent of r), however, it is a function of P only. 

The equation of transfer is of great generality. It does not 
require that radiation is the only mode of energy transfer, it just 
refers to the part of energy appearing in the form of radiation. 


§3. RADIATIVK EQUILIBRIUM. 
RADIATION AND TEMPERATURE 

If 4:7T€ (P)dm denote the heat energy liberated in dm, i.e. the sum, 
net gain within dm of heat energy due to convection and conduc- 
tion plus the radiant energy (of all frequencies) coming from sub- 
atomic sources within dm per unit time, the conservation of 
energy is expressed by the equation 

rZ(T = Itt pedv, (17) 

where V is an arbitrary volume, S its surface and n the outward 
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normal on S, According to (5), integrated through the whole 
spectrum, this is equivalent to the equation 

div ^ = 4:np€. (17') 

It is of importance to realize that (17) or (17') refers only to the 
total radiation of all frequencies. Scattering has no influence 
upon the wave-length; the absorbed radiation is, however, 
re-emitted within other parts of the spectrum, implying, in 
general, an exchange of energy between the different parts of the 
spectrum. 

Integration through all directions r of the equation of transfer 
(after multiplication with /^(a^-for^)) yields, according to (11), the 
equation ^ (18) 

the scattered radiation having dropj3cd out as it should do. On 
integrating ( 1 8) through the whole spectrum we obtain, according 
to (17'), the fundamental relation 

r, (P) = -^-jdc. a, (P) (P, r) dv + e{P). . . .( 19) 

If (19) is multiplied by 47Tdm, the left-hand side becomes, ac- 
cording to (8), the total emission of dm, whilst the first term on 
the riglit becomes the part of the total radiation that is absorbed 
by dm. (19) means therefore that the radiation emitted minus the 
radiation absorbed equals the energy liberated. 

Since (18) is a consequence of the equation of transfer only, 
(17') follows, conversely, from (19). (17') and (19) are therefore 
eipiivalent expressions of the conservation of energy. 

In the case € = 0, used as a first approximation in the outer 
layers of a star, we speak of strict radiative equilibrium. Energy 
is transferred by radiation only. It should be noted that e could 
be negative as well as positive in the general case. When, for 
instance, the element loses heat by convection or conduction, 
subatomic energy being missing, e is negative.* 

* (17'), (19) possess, of course, general validity, as they express the energy 
principle. They are, however, of importance only when radiation is the principal 
agent of energy transfer. It is in this sense that (19) is spoken of as the equation 
of radiative equilibrium. 
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Pure absorption. Gray material. We consider the case 0 -^ = 0. 
Furthermore, let us suppose that the absorption coefficient is in- 
dependent of the wave-length, oi^ = (x, (gray material). This case 
is of special importance because of its simplicity. The fundamental 
equations (10), (16') and (19) become, after integration through 
the spectrum. 


1 

a(P) 


dI{P,r) 

ds 


= J(P)-1{P), 


(20) 


j(P) 


4 / 


/(P, r)dco + 


e(P) 

a(P)' 


(20') 


Monochromatic radiatwe equilibrium. In the case of })ure scat- 
tering, = — we know that the wave-length remains un- 

changed. The scattered radiation is fully re -emitted with the 
same freciuency and distributed through the different directions. 
This state of affairs is often called monochromatic radiative 
equilibrium. The e([uations of transfer (JO), (!(>') take, here, the 

1 dT,(r,r) 
a,(P)' (Is 




.(21) 


■/. ( P, r) = J /, ( /^ r')y,(P-, r', r)d(v' (21') 

In this case, we have, of course, according to (IS), 

cliv5^ = 0. (22) 

It is physically plain that the case of purely absorbing gray 
material in strict radiative equilibrium is formally equivalent to 
the case of monochromatic radiative equilibrium and uniform 
scattering, The absorbed radiation is uniformly re- 

emitted in direction, in the same way as the scattered radiation 
is uniformly redistributed in direction. The equations (20), (20'), 
€ = 0, have accordingly the same form as (21), (21'), 1. This 

formal equivalence is of importance, since it allows us to develop 
the same mathematical theory for the two physically different 
cases. 


Local thermodynamical equilibrium. In the case of pure 
scattering, radiation has no relation to the tem])erature. As soon, 
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however, as absorption and emission play a role, we should have 
information how the radiation is related to the temperature of 
the matter. 

Let ( T) be the intensity of the v-radiation of a black 

body at temperature T, 


B 


V 


2hv^jc^ 

^hv!kT _ 


(23) 


where h signifies Planck's constant, k Boltzmann’s constant and 
c the velocity of light. For the total radiation we have 

rBJv=''T\ (24) 

Jo TT 

a being Stefan’s constant. 

When the radiation is enclosed between black walls, it obeys 
Kirchhoff’s law = (25) 

The condition of being enclosed is, of course, not rigorously 
fulfilled in the stars, since there is always a radial net flux of 
radiation. Eddington’s perfect gas star, however, has such a high 
o})acity that the radiation is practically enclosed, the degree of 
accuracy being higher than in laboratory experiments. Astro- 
physics uses the term ‘local thermodynamical equilibrium’ in 
order to cx})ress that the radiation behaves like an enclosed one, 
i.e. that Kirchhoff’s law holds. Even in the outer layers of a star 
this state of affairs has proved to be a useful approximation. 


§4. THE MAIN PROBLI^ni 

In the physics of the outside of a star, the following problem plays 
an essential role. The radiation field, in particular the radiation 
emergent from the surfiice and its distribution in direction and 
wave-length, and the temperature distribution are to be deter- 
mined when the coefficients of absorption and of scattering as 
well as the law of scattering and e are given. 

This problem naturally represents only a part of the general 
problem of the determination of the whole physical state. The 
quantities regarded as given in the above problem are not strictly 
known in the stars. Actually they enter other more or less known 
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physical laws (containing of course the conditions of mechanical 
equilibrium) which, together with the laws of radiation, should 
enable us to determine the total physical state. The excessive 
difficulties, however, connected with this problem as it stands 
compel us to select partial problems, of which the one formulated 
above is one of the most important. Comparison with observa- 
tional data of the solution of this problem, taken for different 
choices of the given quantities, has led to valuable insight into 
the structure of stellar atmospheres. The following compara- 
tively simple special cases have hitherto found the chief attention 
of astrophysicists. Firstly, the Schuster-Schwarzschild model 
of a purely scattering atmosphere. Schwarzschild was the first 
to approach this model with rigorous mathematical methods. 
Secondly, Schwarzschild’s (formally equivalent) model of a 
])urely absorbing gray atmosphere (€ = 0). The im])ortant limit 
case of infinite optical depth has (by approximate methods) 
been extensively treated by Milne, who also studied the spectral 
distribution of the emergent light. The explanation of the law of 
darkening on the sun’s disk is one of the main successes due to the 
study of those models of radiative equilibrium. Schwarzschild’s 
investigations have, furthermore, led to the result that the origin 
of the Fraunhofer lines in the solar spectrum is due to scattering 
rather than absorption. 

So far as concerns the general solution of the above problem, 
the following remarks are of importance. The equation of transfer 
( IG) can be considered as a first order differential equation for the 
intensity. On integrating it along a ray, with respect to the 
boundary conditions, we obtain the intensity expressed in terms 
of the Ergiebigkeit. If this is inserted into (16'), a linear integral 
equation (due to Schwarzschild) is obtained for the determina- 
tion of the Ergiebigkeit. In the case of pure scattering, a^, = 7 /^ = 0, 
this equation completely determines the Ergiebigkeit. In the 
photospheric layers of a star, however, absor])tion processes are 
predominant, and the unknown emission 'q^ — oL^B^ enters. The 
solution of Schwarzschild’s integral equation determines, in this 
case, the Ergiebigkeit, and therefore the intensity as a ‘ func- 
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tionaF of B^, for every frequency v.* Since = is a 

definite function of B, the intensity appears thus a definite 
functional of B(P), i.e. of the temperature distribution. If 
the intensity thus determined is inserted in the equation of radia- 
tive equilibrium (19), the right-hand side becomes a definite 
functional of B (P), The left-hand side 

r)(P)=roiAP)B,{B(P)}dv, 

Jo 

however, represents a definite function of the point P and of 
B = B(P). We thus arrive at a definite functional — more pre- 
cisely, at a non-linear integral equation for the determination of 
B(P), i.e. of the distribution of temperature. Once B{P) is 
known, we arc able to find and therefore for every frequency. 

It is difficult to solve the integral equation in its full generality, 
and the limitation to simpler special cases seems necessary. If 
the material is gray, i.e. if is independent of v, the equation 
of radiative equilibrium simplifies considerably. If, however, 
varies with v, the integral equation still remains non-linear. In 
gray material without scattering we have the simplest case, since 
the integral equation then becomes linear, and the determination 
of B{P) appears entirely separated from the determination of 
the spectral distribution, the latter being a matter of direct 
integration. It should be noted, that the mathematically more 
general case, where and o-^ are both independent of v, leads to 
the same formal integral equation (linear) as in case = 0. But 
the spectral distribution of light must then be determined by the 
solution of Schwarzschild’s integral equation. 

The greater part of this tract will be devoted to the mathe- 
matical problems arising from the two })articular cases, pure 
scattering, and pure absorption in gray material. The first case 
plays, as mentioned before, a chief role in the outermost layers 
of a star, whilst the second case has found chief application to the 
photospheric layers. 

* *1 his (letermination of is possible by means of the Neumann series for 
the solution of the integral equation. Convergence and uniqueness of that 
solution are treated in § 34. 
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So far as concerns the boundary conditions, the radiation in- 
cident at the surface of an isolated star is zero, while the radiation 
at great depth becomes isotropic. In the case of a close binary 
star, the radiation of the other component must be taken account 
of. Or, in the case of a ‘ planetary atmosphere ’ in radiative equili- 
brium, the solar radiation must be considered. In all these cases 
the radiation incident at the surface is given. Another type of 
boundary condition occurs in Milne’s spherical model of a 
‘ planetary nebula’ in radiative equilibrium. The inner ftice of the 
shell receives not only the radiation of the central star but also 
the light coming from other parts of the inner face. The boundary 
condition appears here as a relation between incident and 
emergent radiation. 


§5. SCHVVAUZSCHILD’S MODELS. MATERIAL 
STRATIFIED ]N PARALJ.EL PLANES 

In the atmospheres of celestial bodies the material may be con- 
sidered stratified in concentric spheres. The radiation field having 
spherical symmetry too, the intensity 1 (F, r) becomes a function 
of the distance a from the centre and of the angle 6 between the 
radius vector and the direction r only, 1 == I (aj6),OS0 In this 

case we obtain 

-j =cohO ^ . 

(Is d(i a vU 

When, instead of a, 0, the variables 

^ = acosO, rj — a^inO 

are introduced, we have more sim])ly 

(U dl 

Neglection of curvature. Let us consider a layer of 

the star. If, the depth being kept fixed, the radius tends 

to infinity, we formally obtain the limit ease of a plane slab. 
Another way of getting this limit case is the following. We set 

^1 ® / , \ ^0 /•/ \ 

* = . P («.. + „ /(^)> 

(l^ (Iq W-Q 

Xq being a fixed quantity, while f(x) is a fixed function of x, 
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0<x<Xq. On proceeding to the limit , i.e. on making the 

shell thinner and thinner and making p -j- cr^) larger and larger 
in the same ratio (at homologous points), we obtain, as seen from 
(16) and from the above expression of dljds, again the limit case 
of a plane slab. This time the radius has been kept constant. It is 
customary to neglect the curvature in the outer layers of a star, 
the material being thus stratified in parallel planes. Let x gener- 
ally denote the depth of the ])oint P below a fixed layer, x being 
reckoned negative for points alcove that plane. The direction r is 
characterized by two angles 0, (f) (we retain the second angle in 
view of later applications), 6 being the angle between r and the 
direction of negative x (outward normal of the slab), 

<j> is the azimuth, i.e. the angle between the plane, containing r 
and the .r-axis, and a fixed plane through that axis, 0 ^ g 2tt. 
We shall only consider the case where the radiation field is the 
same in all points witli the same x, becoming then a function 
of (.r, r) alone. We have 

d8= — sec0dx% (26) 

thus yielding — cos 0 . 

ds dx 


The SchwarzscMhl- Milne rnodel. This is the (jlassical case of 
purely absorbing and gray material in local thermodynamical 
ecpiilibrium. Tt is convenient to introduce, instead of x, the 
optical depth 

T= oipdx 

J —00 

below the surfoce. In the case of strict radiative equilibrium, the 
fundamental equations (20), (20') take then, according to (25) 
and (26), the remarkably simple form 




(27) 

and j*/(T, r)(/aj, (28) 

where, in the case of local thermodynamical equilibrium, 


a Jo 
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If convection, conduction or subatomic energy sources are not 
negligible, we must add the term c/a on the right-hand side of (2S). 
An important fact is that the absorption coefficient appears only 
indirectly, namely through the mediation of the now independent 
variable r. It is to be remarked that (27), (28) determine the total 
radiation field, in particular the direction distribution of the 
emergent light, without any hypothesis about temperature and 
its relation to the radiation. VVe need it, however, as soon as the 
spectral distribution, or the temperature distribution, is wanted. 

Let us introduce the notation 




for the integrals over the hemispheres of all outw ard and all 
inward directions. We set 


TtF^ ~ ^ ^ ^ Gdo), 7tF- ^ ^ ^ ^ 

( 20 ) 

The equations ^ ^ 

I cos ddco= \ I cos d\da> = tt, 


following from , • n in 

° dco — su\0d0d(l>, 

show that F,^_ and F_ rc})resent averages of the intensity within 
the respective hemispheres. Furthermore, we set 

7tF — 7rF^ — TrF_=jI COi^0da>. (* 10 ) 


The flux vector ^ is, of course, normal to the layers, and we have, 
according to (17), e = 0, 

= 77 /<" = const. (.‘11) 


Another important relation, due to Eddington, is obtained if (27) 
is multiplied by cos 0 and integrated through all directions. 
According to (30) and (31) we get 




iv = - / cos‘^ 0dco = Ft -f- const. 


(32) 


The physical meaning of K is that it equals c/tt times the radiation 
pressure normal to the layer. 
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In the absence of strict radiative equilibrium, e ^ 0 , the funda- 
mental equations are 0 ^ 


cos 0 ^— = 1 — J, 

CT 


.(33) 


J = -’ (ldw + -. 

477 J a 

Instead of (31) we have ^ 

dr a’ 

and (32) must be replaced by 


dli 

dr 


= F. 


(34) 

(35) 


(36) 


Another simplification of the model is effected by assuming the 
slab to have infinite oj)tical depth, 

O^T<00. 

This corresponds well to the conditions within a star (large 
opacity in the interior). We formulate then the first standard 
j)roblem, treated approximately by Milne: 

Problem I. I (r, r) ^ 0 and B (r) ^0, r < oo are to be deter- 
mined from (27) and (28), the radiation incident at the surface, 

/((), r)^0; 0>77-/2, 

being given, as well as the flux constant F. 

The incident radiation is zero for ordinary stars; positive, 
however, for close binaries (or in Milne’s model of a planetary 
atmosphere subject to insolation). For an ordinary star, the flux 
constant is derived from the observed effective temperature 2 \ . 
At the surface we have F = F^ , this being an average of the in- 
tensity of the emergent radiation. If this were due to a black 
body, the intensity would be aT^jir, T being the surface tem- 
perature. In general, 

F= (37) 

IT 


defines the effective temperature . 

The relations (27) and (28) being linear in /, J = JS, the solution 
of Problem I will be the sum of the solutions of two partial 
problems: 
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Problem La, Problem I in the case 

/(0,r) = 0, ^>77/2, (38) 

but for an arbitrary F ^ 0. 

Problem T^. Problem 1 in the case 

F = 0, (39) 


In lack of strict radiative equilibrium, we have an 

analogous problem. We confine ourselves to the case (38): 

Problem 11. / (t, r) ^ 0 and are to be determined 

from (33) and (34), the incident radiation being zero, the net dux 
Fq at the surface being given, and e/a being given as a function 
of T. 


The Schuster- Schwarzschild model. This is the case of mono- 
chromatic radiative equilibrium, with e = 0 . The affix v can be 
omitted, the frequency being always the same. On introducing 
an analogous optical depth 


T ~ 



the fundamental equations (21), (21') take, according to (2(3), 
the form 


cos = I (t, r) - J (t, r), 

(jT 


•J (r,!-) =J^j/(r,r')y(r;r',r)dco' (41) 


According to (22), the monochromatic net flux is constant. It is 
to be noted that Eddington’s relation (32) holds too. For multi- 
plication of (40) by cos^ and integration with respect to r yields 

=/<’-- fj (r, /•) COS ffdoj. 
dr IT J 

If, here, the integral is evaluated by means of (41), the integrand 
of the thus obtained r' -integral is seen to contain the factor 

J cos 9y (r; r, r') doi — 0. 

The vanishing is a consequence of (11")> i e. of the fact that y is 
an even function on the r-sphere, cos 6 being an odd function. 


11 M P 


2 
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Problem III. Case of infinite optical depth. I (r, r) ^ 0 and 
J (r, r)^0 are to be determined from (40) and (41), the incident 
radiation being zero and the law of scattering y as well as the flux 
constant F being given. 

Problem IV. Case of finite optical depth, 0 < r < t*. / and J 
are, for a given law of scattering, to be determined, when the 
radiation incident at the outer face is zero and when the radiation 
incident at the inner face, 

/(t*, 0), 0<7r/2, 

is given. 

In the case of uniform scattering, y = 1, Prol)lem II I is formally 
the same as Problem la. The case y=l and /(t^, 0) = const., 
0<7r/2 (black body as a background), is the classical ])roblem 
treated by Schwarzschild in a fundamental memoir. Problems 
III and IV have hitherto not been discussed in full generality. 
The existence and uniqueness of the solution will be discussed in 
Cha])ter iv. Fiirthermore, limits for the solution will be given 
that are independent of the law of scattering. They are the same 
limits which were found by Schwarzschild in the case of uniform 
scattering. In this classical case, considerably narrower limits 
will moreover be found. 

§G. FIRST CONSPJQUENCES. THE RADIATION 
FROM GREAT DEPTH 

We first draw some simple conclusions from the positivity of the 
intensity and of the Ergiebigkeit, confining ourselves to the two 
standard cases formulated in the preceding section. We suppose 

^ 0 for all sufficiently large r. 

(35) shows then that F (r) never increases, for t sufficiently large. 
From (36), we have 

K(T)-K(0)=j^F(t)dl. (42) 

Since, according to / ^ 0, if is never negative, (42) implies 

j^Fdt^ -F(0). 


(43) 
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This proves that jP ^ 0 holds for all large enough r, i.e. that the 
net flux is outward for these r, for the opposite assumption would, 
according to the decreasing of F, imply that the left-hand side 
of (43) tends to — cx) as t->oo. These considerations show also that 


the limit 


Fao = lim -P (t) ^ 0 


exists. From this fact and from (42) we obtain the important 

relation I ^ 

Foo = Iim— cos^ 0 1 (r,r)dio, (44) 

T — QO J 

The generality of this equation must be emphasized. In strict 
radiative equilibrium, e = 0, F^^ is of course the flux constant jP. 
(44) holds in the Schuster-Schwarzschild model too, because 
Eddington’s relation was found to hold in that case. 

From (35) we now get 

7rF^ = 7rPoo -f Itt f -cZt;Foo^0. (45) 

J T ^ 

The second term, taken for r = 0, 

47rJ^lrfT = 47rJ^ pedx, 

represents, as it should do, the energy liberated within a normal 
cylinder through the whole slab, having the unit area as its 
cross-section. As a special case of (45), we note that 

Foo = Fo~4r"-dT. (45') 

Jo ^ 

It is thus seen that Problem II additively decomposes into two 
partial problems, the first one being Problem II in the case 
JPoo = 0, while the second one is simply Problem la, F=Foo. 
We therefore confine our attention to 


Problem II a. Problem II in the case 

Jo a 

i.e. in the case where the net flux at the surface equals the energy 
liberated per unit time within a cylinder of cross-section one. 
The solution of the general Problem II is obtained in adding 
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to the solution of Problem II a the solution of Problem la with 
an arbitrary F — Frj^ ^ 0. 

On solving the equation of transfer (40) with respect to the 
intensity, we get 

e'’’'“c9/(7.',r) = c J {t,r)dL 

(46) 

It may seem surprising that no boundary condition for t->oo 
appears in the formulation of Problems I, II and III. The su])- 
position, however, that / and J be non-negative quantities, 
allows us to dis})ense with it. 

Theorem I. We have 

I (t, r) = sec e j " ,/ {t, r) dt; 0 g 0 < 7 r/ 2 , 

( 47 ) 

i.e. the radiation coming from the deeper layers (0<7r/2) is 
solely due to the Ergiebigkcit of the material within the slab. 

Proof. Considering in (46) the case 0<7r/2, t<t', we infer, 
according to J ^ 0, that, for any fixed direction r, 6< 7t/2, 

I {T\r)\ 0g(?<7r/2, 

decreases w ith increasing r'. According to / ^ 0 this implies the 
existence of the limit 

/(y')=lim e ^ 0, 

t'—co 

for We are thus allowed to i>roceed to the limit t'->oo 

in (4G), 

^00 

/ (r, r) - i (r) + sec 0 J j i^t^y)dcDj 

for 0 < 7r/2. "Ilie first (non-negative) term on the right-hand side 
represents the radiation not due to the Ergiebigkeit of the slab. 
Such radiation cannot, however, exist in the steady state. We 
have from (48), according to sec0^ 1, 

I (t, r) ^ i (r) ^ 0, 


thus yielding 
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This is, according to the finiteiiess of Fao, compatible with (44) 
only if ^ 

J icos-0rfco = O, 

i.e. if i (r) vanishes for almost all directions /* (in the sense of the 
Lebesgue measure). Values different from zero in a set of measure 
zero are, however, physically meaningless. It should be noted 
that we cannot measure intensities, but only amounts of energy. 
Such an amount is always represented by an integral J / do) ex- 
tended through a finite cone of directions. The values of / in a set 
of directions of measure zero are here, of course, without inttuence. 


§7. THE TXTE(;HAI. EQUATIONS OF THE PROBLEMS 

The reduction to integral e(| nations of boundary value problems 
of the theory of radiation goes back to Hilbert and Schwarzschild. 

Mihie^s integral equations of Problems /, II. On setting t = 0 
in (46) and writing r instead of r' afterwards we find, for 0 > 7r/2, 
the radiation coming from the upper layers, 

i(T,r) = e-^l*‘’"‘^^l/(0, 0+ |scc(?| f J(/)rfL 

Jo 



When (48) and (47), J = J(t), arc inserted into (34), a linear 
integral equation for J (r) is obtained. In order to find its kernel 
we introduce, according to Schwarzschild, the functions 




■" 9dio 


.(4!)) 


= 2 ^/ 


e-r|scc0i|ycc2 


On setting s= | sec0 | they appear, according to doj = i^inOddd(l}, 
in the well-known form 


{x) = J ds. 


.(50) 


Insertion of (47) and (48) into (34) yields two essential terms, the 
first one being i ^ 

a 47rJ _ 
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due to the energy liberated and to the incident radiation, while 
the second term is the sum of two double integrals 


1 

4-77 



sec 0 6“" ^ J {t) dt doi 



sec 0 I I I J (t) dtdoj, 

(51) 


Tlie integrand being measurable and non-negative, we may inter- 
change the order of integration, obtaining thus, according to 
(49), 71=1, , 

A(./),s-J^ J{t)E{\r-l\)dl, (52) 

for (51). The use of a brief symbol A for the linear integral operator 
is convenient for the treatment of the problem. We thus obtain 
Milne’s first integral equation for the determination of the 
Ergiebigkeit J (r) = B (r), 

J(t) = A(J)^ + ' ■ + /- I* ...(53) 

a 477 J _ 

A formal disadvantage of (53) is that it does not contain the net 
flux TrF{r). An equation containing F is, according to Milne, 
obtained by inserting (47) and (48) directly into (30). According 
to (49), 71 = 2, we thus get Milne’s second integral equation 

ir(T)=^J cosfle”’‘'®®‘^^>/(0,r)Qfct> 

+ 2 J J(t) /?2 (t. - r) dt J (0 ("r - 0 dt, 

(54) 

containing all the data given in Problems I, II. (54) goes, of 
course, together with (45). 

If (47), (48) are inserted into the integral (32) representing K, 
we obtain, according to (49), 7i = 3, the equation 
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Since the first term on the right-hand side is not greater than 

A"q = - J cos‘^ dl (0, r) doj, 

it follows that (44) can also be written in the form 



For the actual solution of Problems T, II, it is more convenient to 
start from the integral equation (53) together with (56) and (45'). 

Theorem II. Any non-negative solution B(t) of (53), that 
satisfies (56) where F^o i« determined by (45'), yields together 
with (47), (48) a solution of Problem II, with given energy 
liberated, given incident radiation and a given surface flux 
ttFq ^ 4:7t jpedx (integrated through the slab). 

Corollary 1. In case of Problem la we must, in (53), set 6 = 0 
and I (0, r) = ^ 

Corollary 2. In case of Problem lb we must set 6 = 0 in (53) 
and F = = 0 in (56). 

Corollary 3. In (;ase of Problem Ur we must set I (0, r) = 0, 
6 > 77 / 2 , and Foo = 0 in (56). 

The proof of the theorem is plain from the above. 

The integral equation for Problem HI. Since the incident 
radiation is supposed to be zero we have from (46) 

/ (r, r) = I sec 0 1 J’’ c-Mt-O iHocO I J ^,5. e> IT 1 2 

(•'57) 

for the radiation coming from the upper layers. From (IT), (11") 
and from (41) we infer that the Ergiebigkeit is always an even 
function of direction, J (r, -r)^J(T,r). Let us introduce the 
linear integral operator 

A J r')J(t, r')dtdoj', 

(•'iS) 

where H is defined by 

H (T,r;t,r')= | sec0' j (59) 
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for any two directions r, r\ is thus the kernel of the operator 

A. Tn the simplest case of uniform scattering, y = 1 , this operator 
reduces to the simpler Milne-operator, J being independent of 
direction. For an arbitrary law of scattering, however, no 
simplification is possible. The meaning of the general operator is, 
that, on inserting (47) and (57) into (41), the linear integral 

.7(T,r) = A(./),., (60) 

is obtained for the determination of the Krgiebigkcit. In analogy 
to Theorem II we have in the ])rescnt case 

Theorem [11. Any non-negative solution of ((Jl) satisfying 
(44), F,j^ = F, yields together with (47) and (57) a solution of 
EVoblem 111. 

The. integral equation for Problem TV. Since the incident radia- 
tion vanishes, (57) gives the intensity for 6 > 77 / 2 . The intensity 
of the radiation from the deeper layers is, however, given by 

= (T^r)-fsec^ J {t,r)dt; 


<9 < 77/2 (61) 

We define, in this case, the integral operator 

L{J)r.r m 


by replacing, in the integral in (58), the upper limit 00 of integra- 
tion by T*. On inserting (57) and (61) into (41), we obtain the 
linear integral equation 

J (r, r) = L (,/)^,, + J (r*, ?•') dcj' 

■ ' (63) 

for the determination of the Ergiebigkeit, 0 ^ r ^ r^. 

Schwa rzschild’s integral equation for uniform scattering. Since, 
in this case, J is independent of direction, the operator L{J)^ ^ 
simplifies to ^ 

i=/v(./), = ^,J^ J{t)E,(\T-t\)dt, (64) 

and the integral equation (63) becomes 

J{t) = L ( J)^ + J / (t*, 6') do}'. 


( 64 ') 
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The classical case treated by Schwarzschild is 

/ (t*, r) = /* = const. ; 0<7r/2. (65) 

In this case the second term in (64') becomes, according to (49), 
n = 2, simply ^ (t* - t). (66) 

For later purposes we add the remark that in the case 

/(r*, r) = cos<?; 0<7r/2 (67) 

the second term in (64') becomes, according to (49). n = 3. 

\Kdr*-r). (6S) 


Positivity of the kernel. An obvious as well as imjjortant pro- 
perty of all the integral operators introduced above is that their 
kernels take only positive values. We therefore call A. L positive 
operators. This positivity occurs in all boundary value problems 
of the theory of radiation. A?! evident consequence of that 
property is the 


Lemma. For any functi<^n <I) (/, r') ^ 0 we have A (d>)T,r > 0 
everywhere, unless <!> vanishes identically. 'I’ho vanishing is, 
strictly, to be understood in the sense of the theory of the Lebesguo 
measure, i.e. vanishing up to an inessential (r, r)-set of measure 
zero. 

We may express the lemma also in the following form. <[> ^ 0 

and 0^0 implies A(<I))^,.>0 for all t, r. Another equivalent 

formulation is that ^ ^ . 

T r ^<P 


implies the strict inequality 

A('F),,,>A(0),,, (69) 

for all T, r. The equivalence follows from the linearity of the 
operator, in particular from 

A('F) = A((&) + A('F-(I>). 

The same thing is, of course, true for the operator L. 
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§8. SOME PROPERTIES OF THE FUNCTIONS E^{x) 

For the discussion of the simpler integral equations we need some 
simple properties of the functions x‘^0, n^O. VVe have 

= while Ei{x) becomes logarithmically infinite as 

For /«.> 1 all functions (*r) arc, however, continuous at 
x~(). By partial integration in (50), the well-known recursion 
formula ^ = e--- .f/'J,, (x) (70) 

is obtained. We have, furthermore, 

(71) 

From (50) the ine(|uality 

y<;„,,(.f)<A;,(.r) (72) 

is obtained. Let us compare the right-hand side of (70) with the 
right-hand side of the equation obtained from (70) by replacing 
n with n— 1. We then find, according to (72), 

in-l)EJx)^nE,^,(x) (73) 

for 1, the equality sign holding only at .r = 0. We note that 



If the left-hand side of (70) is combined once with (72), another 
time with (73), we find 

p X p — X 

— T (75) 

x + 71 x+n—l ' 

for ri^l. On applying Schwartz’s inequality to 

(*00/ s n~~l \ / s n-flX 

\ ^ ^ ^ )v 

another inequality 

E^^(x)<E„_i{x) 

tl-hl (^) (76) 

is obtained. From (71) and (76) we get 

d (E„^A 
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■<0; x>a>0 


.(78) 


For later purposes the formula 

A 

dx (:r) 

will be useful. For the proof we note that the left-hand side 
equals, according to (71), 


E,,ix-a) 

E,(x) \ E,(x) E,Xx-a) /* 

By means of (77), the quantity within the parenthesis is seen to 
be negative. 

Some simple integral formulae will, furthermore, be needed. 
From (71) and (74), we have 

£ E, (r-t)dl = y^- (r), JJ* E„ (< - t) = - E„_,, (r* - r). 

^ " ( 79 ) 

Partial integration yields, according to (71), 


j: 


(80) 


lE,^ (0 ^ j - ^,,+2 {t)-tE„ , 1 (t). . 

On combining (79) and (80) we get 

[ t E„ (r -<)(//= f" (r - 1) E„ (<) dt=E,.^{T)- ]- + y 

Jo Jo n-^ i ti 

(«i) 


and 




{t)dt 


= re +T + re ~ 

( 82 ) 



CHAPTER II 

SOLUTION OP PROBLEMS I AND II 

§0. TKM SOI.rTION OF lUiOBLEM Ta 

Millie's integral equation of the Problem la is, according to (38) 
and (53), € = () and J = /I, 

= (83) 

Physical reasoning would lead us to the conjecture that, in strict 
radiative e([uilibrium, the radiation becomes nearly isotropic in 
great depth, i.e. according to (28), that / (r, r)!B (r) I as t->oo. 
(44) would then show that B(t) is asymjitotically linear, 

i = I cos- ddo ) . hm , 

i.e. that B(r) = lFT for large r. We have, indeed, the 

Tukorem IV. Problem la has a solution 

,/(T)=B(T) = fP/(r), (84) 

where / satisfies Milne’s integral equation (83) and the inequalities 

/(T) = T + g'(T); i<5'(T)<l. (85) 

Proof. Prom the definition (52) of the operator A and from 
(79), T* = oo, we have i = A(l),+ U\( t). (86) 

Furthermore, from (81) and (82), n= 1, t* = oo, 

r = A(<),-U’3(r). (87) 

On combining (86) and (87) we get 

r + c = A{t + c),+ l{cE, (r) - E, (r)} (88) 

The smallest value of c that makes the second term ^ 0 is, ac- 
conUng to E^> E^ and E.^ for large t, c= 1. We thus get, 
for all r, y ^ ^ (ggj 

On the other hand, r. =; i is, according to the greatest 

value of c making the second term ^ 0 for all r. We thus have 
A (t) < A (/i)/, fl{r)=:T+^ (90) 


for all T > 0. 
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We now set successively 

A+i(t) = A (/„),. (91) 

The inequality (90) has then the simple form 

/iW</2('r)> (*>2) 


According to the positivity of the operator A, wo have from (92) 
A(/i).<A(/2),, t^O, 

or, with regard to (91), /2 (t) < /., (r), r ^ 0. Successive application 
of the operation A generally shows that 

/.W<A+i(t), t^O, (93) 

i.e. each function of the sequence/,^ lies above the preceding one. 

From /i(t)</(t), t^O, 

we now have, again according to the positivity of A, 
/2W=A(/,),<A(/)„ t^o. 

From (89) wo get therefore / 2 (t)</(t). If the operation A is 
applied a second time, we get 

/3(t)=A(/2),<A(/)„ T.^0, 

whence, according to (^S9), / 3 (t)</(t) follows. In the same way 
we obtain the general inequality 

/„(t)<A(/)^</(t), t^O. (94) 

(93) and (94) evidently show that the limit function 

/(T) = lim/„(T) 

n— 

exists and that it lies between the limits J and /=t+ 1, 

T^O. It follows from well-known theorems that, in (91), we can 
proceed to the limit = oo under the integral sign, /(r) is therefore 
a solution of/=A(/) with the property (85). The smoothing 
property of an integral operator shows, finally, that /(r) is a 
continuous function. 

On setting B (r) = cf (r), 

it remains to find the relation between c and the flux constant F 
From/=T + g and from (81), (82), n = 3, t* = oo, we find 

f fif)E^{\r-t\)dt=lT>rE^{T)-k-[ q{t)E3{\r-t\)dt. 

Jo Jo 
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The last term being, according to 0<^<1 and to (79), 7^ = 3, 
T^=:oo, less than §, we infer from (56), J = cf, = that F 
equals fc, i.e. that (84) holds, q.e.d. 

From /=A(/) and from (87) we infer that the remainder 
function q (r) satisfies the integral equation 

q{r)^\{qME.,{T), (95) 

Let us, now, insert (84), f=T-{-q, into the second integral 
equation (54), where e and the incident radiation are to be omitted. 
Making use of (81), (82), t^ = oo, we find that ^(t) also 

satisfies the integral equation 

/ {t) (t -r)dt- JJ q (/) (t -t)dl = E^ (t). . . .(96) 

On inserting (84) finally into (55), we o])tain by means of similar 
computations, and according to (32), the equation 

f (/ (^) (I T — ^ I ) A5 (r), (^"^) 

J 0 


a being a constant, the determination of which will be post- 
poned. 

From (22), (34) and (87) we find the temperature T as function 
of the optical depth r, 



B{r) 

y 


= U{ry 


(98) 


Formulae (47), (57) for the intensity can be written in the form 


e^'^J {T + 8Cosd,r)ds; 0<7r/2, 


.(99) 


rr 1 sec0 1 

^(^>0= e~^J {T + scos6yr)ds; 0>7r/2. ...(99') 

J 0 

In the present case we have I == I (t, 9) and J (t)== B{r). On 
inserting (84), (85) into (99) we get 


I -IP T -h cos 0 + J g (r -h 5 cos 0) ; 


0 ^ 0 < IT 12, 


.( 100 ) 
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for the radiation coining from the deeper layers, in particular, 


I (0,0)=:^ IF 




(101) 


for the radiation emerging at the surface. 

It is to be emphasized that the emergent radiation is in- 
dependent of the particular behaviour of a within the slab, 
provided that the material is gray. Since Problem I II for uniform 
scattering formally coincides with Problem la, we may also say 
that the emergent radiation (for a given frequency) does not 
depend on the course of It should be repeated that (101) 

also holds independent of any assumption about tem23erature. 
In the case of the sun, (101) rei^rescnts the law of darkening on 
the solar disk, in good agreement with the observations. 


§10. ok TliK SOLUTION 


The uniqueness of the solution of Problem I a follows from the 
Theorem V. A non-negative solution of the integral equation 
(83) is necessarily of the form B (r) = const./ (r). 


Proof. Let B(t) be a non-negative solution of (83). We set 

6 = greatest lower bound of (Ib2) 

JKV 


and B* (r) = B{t) — bf (r), 

/(t) being the solution found in the jneceding section, B^ is then 
also a solution of (83), with the j)roperties 

J^*(t)>0, (103) 

Gr.l.b.of (104) 

j(v 


All we have to prove is that B* (r) identically vanishes. 

Prom (104) and from the definition of the greatest lower bound 
we infer the existence of an infinite sequence of numbers r^, 
v= 1, 2, 3, ..., such that 
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Here, we must distinguish between two different cases, (a) the 
Ty, have a finite limit point t^, (6) we have t^->oo as v->co. 

Let us first treat case (a). It is always possible to pick out of 
the a subsequence converging to r*. No generality is lost in 
assuming that as v->co. From (105) and from /( t ^)>0 

lim7i*(T,) = 0, 

V~r/j 

whence, according to the integral equation, 

lim A(^*)^^ = 0 

l/-=CO 


follows. As a fumition of t, A(B*)^ has the property of lower 
scmicontinuitv, the r)roof of which fact may be T)ostponed, 
IimA(/^X^A(/i^),/. 

V--CO 

This shows with rcs])ect to (103) and according to the positivity 
of A, that B* (r) ~ 0. The lower semicontinuity is proved in the 
following way. According to (103), we have, for 

r* - 8 < T < T* 4- 8, 

A ( /i), a I f ‘ \ (0 E, (r -l)dl + \r B (t) E, (I - t) dt, 

‘“Jo 

(106) 

8 being a positive number less than r^. (In case t* = () the first 
term is to be omitted.) The left-hand side is, now, a continuous 
function within (r^ — 8, t'*^ 4-8), the logarithmic singularity of 
being excluded. We have therefore 

liin A ( B)^ ^ J r * /}/!/, (r* - 1) dt + J r ” BE^ (t - T*) dt. 

(107) 

Since tliis inequality holds with an arbitrary 8, we find for 8 -> 0, 
ImiA(B)^^ A(«)^., 

T->T* 

i.e. the lower semicontinuity. 

Case (6). Omitting the star in B* we have, according to 
/(t) -t, to prove that a solution B{t) of (83) satisfying 



T, 


(108) 
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vanishes identically. According to B^O and to E.^< we 
have /% cx) 

I BE^{\r-t\)dt<2\(B)^^2B{T) (109) 

B{t), being a non-negative solution of (83), determines a 
solution of Problem la. The corresponding net Hux F is not 
known; however, we know from §(>, that necessarily i^^O. 
From (56), J ~B, and from (108), (109) we infer that 

thus yielding F = 0. On inserting this value into (54), J—B, 
incident radiation = zero, wc find for r = 0 




this being compatible with B ^ 0 only if B (r) ~ 0, which com- 
pletes the pro(jf of the uniqueness. 

Corollary. Every solution of (83) with a finite lower bound 
has the form const, /(r). 

Proof. In the solution f3 (r) -h rf (r) the constant c can be chosen 
so large that this solution is everywhere positive. Theorem V 
then completes the proof. 

§11. A DIFKERMNTIATION FORMULA 
Foi^the continuation of the discusvsion of Problems I and II, a 
study of the more general operator 


( 110 ) 


and of the related integral equation will be useful. We introduce 
the notation /• x. 




g(t)h{t)dt. 


The kernel of (110) is evidently symmetrical. The symmetry can 
also be expressed by the general relation 


(g,A{h)) = {li,A{g)) 


n rc 

0 


li {\T — t\ )(j{t)h{T) dtdr. 


H M P 


3 
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We start with an important property of the particular operator 
(110). For this purpose let us suppose, that H (x) is continuous 
for j:; > 0 and satisfies the conditions 

11 = x->0; H {x) = O x->oo. 

( 112 ) 

Let h (x) be continuous for a* ^ 0 and continuously differentiable 
for X > 0. h (x) is, furthermore, supposed to have the properties 

^ A'(:r) = 0^1og^ j, x~^0; h' (x) 0 x->co, 

(113) 

with a suitable 0 < 1 . Under these conditions, the formula 

^^A(/0, = A(A'), + A(0)//(t) (114) 

holds. 


Proof. Wo have 


.(hl=r. r JI(t)h(r-t)(U+ {''n(f)h(r-^t)clt (115) 

Jo jo 


On writing Si (r) for the first integral, we find 

(t + - S, ^ I r ' + 8 _ 


-) I r’"'* 

^ = H(f)h(r + 8-t)dt 


.(US') 


Here, the first term obviously tends to h (0) 11 (r) as S -> 0 . In the 
second term, however, limit process and integration may be 
interchanged. The proof of this follows the same line as the proof 
that a Newtonian potential can be differentiated under the 
integral sign, even when the moving point lies within the mass. 
That the second integral in (115) may be differentiated under the 
integral sign, follows from a well-known theorem, according to 
which this is always permitted within a certain r-interval if the 
T-derivative of the integrand lies, in absolute value, below a 
fixed integrable function of t. This test applies here since, accord- 
ing to (112) and (113), 

1 // (^) h' (T-\-t)\< const. 
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holds for all sufficiently large t, the constant being independent 
of r. Within an arbitrary finite r-interval, having t = 0 as an 
outside point, the suppositions of the above theorem are fulfilled 
because the absolute value in question is less than const. 
for large t and less than const. | log^ | for small t, 

§12. ASYMPTOTICALLY LINTP^AR SOLUTIONS 
OF/=A if) 

Since the formula of § 1 1 will not be used until § 17 we abandon the 
hypothesis (112) about the kernel. Let us consider the integrals 

H, (x) = II (x), (x) = £” H„ (X) tlx. 

We now suppose that the kernel is positive, 


H(x)>0, (1 10) 

and that H (x) olH, (x) (117) 


holds with suitable a > 0. All these assumptions are fulfilled in 
the case (52), with II — .V/i/x. (117) implies the existence of //g, 
and automatically of all If ,^ . 

The integral equation ^ (/)^ (118) 

will be studied hereafter by means of Fourier integrals. It is, 
however, worth while to consider it first fmm the same, physically 
more natural, point of view as in § 10. (IIS) can have solutions 
of various asymptotic behaviour for r large. Under the particular 

hypothesis // 2 { 0 )=i, (119) 

for instance, we have the 

Theorem VI. Under the suppositions made above about the 
kernel, in particular (119), (118) possesses a solution 

f{r) = T + q{T), 
where a is taken from (117). 

Proof. Ill Milne’s case, H=hEi, (119) is, according to 
Ho (0) = ^2 (0)/2 = J, certainly fulfilled. We first note that 

J (t) dt = xll,,^^ {x) + ^2 (^) (1^0) 

3 2 
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Taking account of (119), we find, after simple calculations, the 

relations 1 = A(l), + ^f 2 (T) (121) 

and r = A(t)^-II^{r). (122) 

According to (117), we obtain from (121) and (122) the in- 
equalities A(r)<A(/.)„ /.-r, 

1 


and 


/(t)^A(/)„ /=t + - 


This shows, in the same way as in § 10, that (118) has a solution 
between t and t+ 1/a, q.e.d. 

From (118) and (122) we find the integral equation 

q{r) = A(ql + II,{T) (123) 

for the remainder function q(r). We must now remember that 
the above solution /(t) of (118) was found to be the limit function 
of a sequence //j (r) defined by the recurrent relations 

On setting = q ^^ , we find 

<7«+i('r) = A(g„)^ + /?3(T), 5-1 = 0, 

thus yielding 3 (t) = S A’' , 

0 

being the ntli iterate of the operator A, A® (/) =/. In other 
words, q (r) is represented by the Neumann series of the integral 
equation (123). 


Definition. If the solution of a linear integral equation is 
represented by the Neumann series, we call it the W-solution of 
that equation. 

The function q (r) of Theorem VI is accordingly the A’-solution 
of (123). 


Lemma 1. A nowhere negative solution of the inhomogeneous 


equation 


^ (t) = A {(j>)^ -h 0 (t), 0^0 


.(124) 


is always the sum of the V-solution of (124) and a nowhere nega- 
tive solution of the homogeneous equation. 
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Proof, From (124) we have, according to the positivity of A, 

(r) = S A*' (O), + A^^+i ((^), > S A*' ((D), . 

0 0 

Proceeding to the limit n = co, we see that the existence of a 
non-negative solution of (124) implies the convergence of the 
iV'-series of (124), this series consisting of nowhere negative terms. 
The lemma follows then from the inequality 

<^(t)^SA-(4)),. 

0 


Lemma 2 . A continuous solution <7 of (118) with the property 
glf^c, T->oo, necessarily equals cf. 

Proof. From the supposition made in the lemma we readily 
infer that the function ^ (r)f(r) 

attains either its maximum or its minimum value for a finite 
value of T, say r^. Without loss of generality we can confine 
ourselves to the case of the maximum value M, In this case we 
have ^ ^ ^ 


The left-hand side is here a solution of (118). According to the 
positivity of A, however, a non-negative solution of (118) can 
only vanish somewhere if it vanishes identically. It follows at 
the same time that M — c, q.e.d. 

Let us now apply these simple results to the equation ( 121 ), 
having 1 as a solution. All requirements of Lemma 1 are fulfilled 
in this case. 1 is therefore nowhere smaller than the A-solution 
of ( 121 ). The A-solution being continuous everywhere, the 
difference is an evidently bounded and continuous solution of the 
homogeneous equation (118). If this solution is identified with the 
g (r) of Lemma 2 , / (r) being the solution of Theorem VI, we find 
from that lemma that g (r) vanishes identically. 1 is therefore 
the V-solution of the equation ( 121 ). We collect these facts in the 


Lemma 3. q (r) of Theorem VI is the V-solution of (123). 1 is 
the iV-solution of ( 121 ). 
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§ 13. AN AUXILIARY THEOREM WITH 
APPLICATION TO PROBLEM II 

Let A, for a moment, be an arbitrary symmetrical operator. For 
the solutions of the simultaneous integral equations 

= <^2 = A(<^2)+«>2, (125) 

we find purely formally 

(^1 > ^ 2 ) “ (^1 5 ^ 2 ) ~ (^ 1 ) ^ (^2))> (^2 > ^ 1 ) (^2 j ^ 1 ) “ (^2 > ^ 

and according to the symmetry of A, 

{<l>M = {cf>M (126) 

This formula is, of course, always correct when the integration 
interval is finite and when the kernel is continuous. The proof 
fails, however, in all of our cases, since here (^ 1 ,^ 2 ) 1® always 
infinite. Nevertheless, formula (126) holds also for certain 
singular integral equations, of course under certain restrictions. 

Auxiliary Theorem. Let A be a positive and symmetrical 
integral oj)erator. Furthermore, let us suppose that 
Og ^ 0 in (125). Formula (126) then holds, if in (125) and b® 
the A-solutions. 

Proof. The iterated operators A*' are well known to be sym- 
metrical once A has this property. From 

00 00 

0 0 

we have 

((Dg , c^,) = S (cpg , A*' (CD,)) = V (o,, A*' (cDg)) = (cD, , <^g). 

0 0 

The termwise integration of an infinite series of functions, applied 
here, is readily seen to be justified, all terms being nowhere 
negative. 

We now return to the special operator (110) and prove the 

Theorem VII. The necessary and sufficient condition that, 
under the suppositions (116), (117), (119), 

(t) = A (<f>)r+^ (t), cD ^ 0 


(127) 
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have a nowhere negative solution <f), is that 


be finite. 



(S>(t)dt 


Corollary. Under the same assumptions, O however having 
any sign, (127) has always a solution if 

be finite. The ^-series converges absolutely. 

Proof. According to (117) and to the fact that H 2 {x) is a 
decreasing function, we have 

H — {t — T)> (i), 

thus yielding, with respect to <^^0, 

(128) 

Furthermore, wo have II — olH^ (t — t) and, since (x) lies 
between positive limits within (0, r), 

V <jy{t)H{T-t)dt'>coxi^t. V <f>{t)H^{t)dt, (129) 

Jo Jo 

the constant being positive and depending continuously upon t. 
A must, now, be finite for some values of t. Hence, and from 
(128), (129), we may infer that (^, Hg) is finite. According to 
Lemma 1 this holds a fortiori, if ^ is replaced by the A-solution 
of ( 127). On applying the auxiliary theorem to the two equations 
(121) and (127) and on taking account of Lemma 3, we find that 

( 1 , 0 ) = (^, H^)<qo. 

The condition is sufficient. On setting 

<^,(t) = SA>'(0)„ 

0 

we infer from the symmetry of A*' and from (121) that 
(<!>, , ^ 2 ) = S (A*' (O), H^) = S (O, A*' {H^)) 
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proving thus that , H^) lies below a finite limit independent 
of n. The fact that (f),^ form an increasing set of functions 


^fi = ^n+l ’ ^ 1, 2, 3, . . . 

implies the existence of the limit function ^(T) = lim^^(T), and 
we infer that (<{>„, II,) 

is a finite quantity. At the same time we readily infer that 
7/2) = ( 0 , 1 ). 

The corollary simply follows from the theorem, since an 
absolutely integrable function O (r) can be written in the form 
O^— Og, Oi and Og being both non-negative and integrable. 

Application to Problem II, Milne’s integral equation (53) is in 
the case of Problem II, J = B, 


= (130) 

a 


Supposing that everywhere e ^ 0, we find from Theorem VII that 
the necessary and sufficient condition that (130) have a non- 
negative solution is, that 

^00 

(131) 

JO a 

i.e. the energy liberated ])er unit time in a normal column of 
cross-section one, be finite. We prove, furthermore, the 

Theorem VIII. Suppose that € ^ 0 everywhere. The necessary 
and sufficient condition that Problem Ha has a solution is that 
(131) be finite. The uniquely determined solution is given by the 
A'-solution of (130). 

Proof, Only the assertion concerning the A-solution of (130) 
needs to be proved. We remember that Problem II a is the special 
case of Problem II where the surface flux ttFq satisfies 

n = (132) 

It is therefore to be proved that the A-solution B{t) of (130) 
automatically satisfies (132). Prom (54), t = 0, J=B, taken 
without incident radiation, we find 

F,^2(B,E,), 



SOLUTION OF PROBLEMS I AND II 


41 


B being the N -solution of ( 1 30), we get from the auxiliary theorem, 
applied to (121) and (130), 

which completes the proof of (132). For any other non-negative 
solution B of (130), however, the left-hand side of (132) would 
be greater than the right-hand side, 

i’o>4 (133) 

Jo a 

In order to obtain the solution of the general Problem II, we must 
remember formula (45'), giving Fco when Fq is given. 

Theorem IX. Suppose that € ^ 0 everywhere. The solution of 
Problem II is uniquely determined by 

5(t)=/7(t) + PV(t), 

Foo being determined from (45'), B being the iV’-solution of (130), 
and/(T) being the function of Theorem IV. 

The uniqueness of the solution follows easily from the fact that 
any non-negative solution of (130) is the sum of the ^-solution 
of (130) and a non-negative solution of the homogeneous equa- 
tion, the latter being characterized by Theorem V. If ttF denotes 
the flux in Problem Fla, we have 

F^=I\ + F^, = F^^O (134) 

If the assumption e^O (everywhere) is abandoned, similar 
theorems hold. But we have then to formulate the general con- 
ditions under which the solution becomes non-negative as 
required in the physical problem. 

§14. OTHER APPLICATIONS. 

THE BOUNDARY TEMPERATURE 

We set (1 35) 

where/(T) is the solution of (118) mentioned in Theorem VI, and 
G'8 (t) = ^ jy II it) f(T + 8- t)dt, (136) 



42 


SOLUTION OF PROBLEMS I AND II 


S being a positive number. From the integral equation (118) and 
from the identity (115), we find 

!78(T) = A(grs).+ Gfs(T). (137) 

Proceeding to the limit S -> 0 we get, purely formally, the integral 

equation for the derivative/' (r). However, as we do not know if 

this derivative exists, wo had better use (137) as it stands. 

In order to prove that is the iV'-soIution of (137), we first 

note that w 

gr8(T) = SA-((7a),+ A-^M9^s)r (138) 

0 

Keeping 8 fixed we infer from Theorem VI that g^ (r) is a bounded 
function of r, 1 | < G, whence 

According to Lemma 3, we have 

a^+M1)=4 s Ay(E^), 

71+1 

showing thus that the remainder in (138) tends to zero as / 2 .->oo. 
( 7 s(t) is therefore the A-solution of (137). 

Theorem X. The value of the solution/(T) of (118), given in 
Theorem VI, is determined by a simple formula, 

/o = V2//7(()]. 

Proof. Applying the auxiliary theorem of §13 to (121) and 
(123), we obtain = = H,(0). (139) 

Since, by partial integration, we find altogether, 

according to /= r + q, ^ 2^, (0). (139') 

Since, furthermore, g^ is the A-solution of (137), the auxiliary 
theorem can be applied to (123) and (137), yielding thus 

{QyG8) = {g8,H^). ..-....(140) 

On the other hand, we find by means of (135) 

(j/s , ^ 3 ) = - 1 J */(<) (0 +/g7(o 

The fraction in the second integral lies between the limits 

H^(t) 
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because = —^ 3 ' is a decreasing function. We may therefore 
proceed to the limit under the integral sign, finding thus 

lim (grs. ^ 3 ) = (*>) + (/, H^) (141) 

8=0 

Furthermore, we have from (136) 

lim(?s(T)=/of? (t). 

8=0 

On inserting the right-hand integral in (136) into the left-hand 
side of (140) it is seen without difficulty that we may proceed to 
the limit S ->-0 under the integral sign on the left-hand side of 


]imiq,Gs)==fo{q,H). (142) 

8=0 

Finally, we have from (123), t = 0 , according to (7o=/o> 

f, = (q, 11) + 11,(0). (143) 

Combination of (140), (141), (142) and (143) leads to 


whence, according to (139'), the required formula follows. 

Lemma 4. Let H (x)IH,{x) be bounded for x large enough. 
Under the condition ( 112 ), the solution /(r) of (118) has then a 
positive derivative /' (r), being the W-solution of the integral 
equation /'(t) = A(/%+/o//(t). (144) 

Proof. First we note that this applies to Milne’s case, H = | , 
since EJE 2 ~^^ x->oo. Since the lemma is only needed for the 
proof of Theorem XI we shall content ourselves with brief in- 
dications of the rather lengthy proof. The main part of the proof 
consists in establishing an inequality for the iV-solution of the 
integral equation 

rj(r) = A{g),+f,H(T). (145) 

The proof is complicated by the possibility of an infinite H ( 0 ) 
(this is just the case with Milne’s equation). It is therefore con- 
venient to split 11 (r) into two parts 

/,^=<D,+a)„ 

J 0<T<o, 

In ^ 


where 


(140) 

(147) 
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<I> 2 (t) is zero for small t and coincides with (r) for large r. 
We consider then the two equations 


= A (0,) , </>, = A (^ 2 ) +^2 (148) 

separately. According to the supposition made in the lemma we 
have ^2 i'^) < const, //g (r) for all r. The iV^-solution of the second 
equation exists therefore as seen by comparing it with the 
equation (121). At the same time we find 

0 < ( f)2 (t) < const. 

for all T. As to the first equation, we begin with considering the 
expression 


A(0 




II (r — t) H (t)dt, T>a. 


We note that A(<1)^) is a continuous function of t. This holds, 
according to (112), also at r = 0. For r large we find, according to 
the hypothesis made, and since H 2 decreases. 


A (Oi)^ < const. J II2 {T-t)H (t) dt < const, //g (r — a), . . . ( 149) 
the constant being independent of a. Now we have, for t large, 
7/2 (t — a) — //g (r) = f II dt 

J T-a 

< const. H 2 dt<a ,cowt,H 2 {T — a), 

J T—a 

Since here the constant is again independent of a, we may choose 
once for all a so small that a . const < J . For this value of a, we 
finclthen (r - a) < 

thus yielding, with respect to (149), 


A (OJ^ < const. H 2 (t) 

for T large. The left-hand side being continuous for all r, we can 
choose here the constant such that the inequality holds for all 
T ^ 0. Comparing again with (121), we find that the JV-series 


SAMO,), 

1 

converges for all r. On adding here the term r = 0, namely , we 
obtain the complete iV’-solution of the first equation (148). 
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Altogether we see that the A’-solution g{r) of (145) satisfies an 
inequality ^ (t) <// (r) + const. 

From (137) and (145) we now find 

(/8 -gr = A igs-g) + Soti)- 

The second term on the right is easily seen to have the following 
properties as 8->0. It is 0 |log for all small t, uniformly with 
respect to 8, while Os(r)-fj,H (r) 

^2(t) 

tends to zero uniformly for all t ^ tq , Tq being an arbitrary number 
greater than zero. A similar splitting leads then to the result that 
g is the limit of ^3 , i.e. that g (t) =/' (r). At the same time we find 

an inequality /' (r) < (r) + const. (150) 

The above considerations show that/' is the iV-solution of (144). 


Theorem XI. Suppose that <I)(t) is continuous at t = 0. 
Under the hypothesis made in Lemma 4, the iV-solution of the 
integral equation (127) has, at t = 0, the value 

( 151 ) 

Jo 

Proof. Since, under the stricter hypothesis about H made in 
Lemma 4, the derivative /' (r) exists, we can dispense with the 
complicated formulae (135), (136), (137) and apply the auxiliary 
theorem of § 13 directly to (127) and (144), getting thus 

(<!>,/') =/o(0,«). 

From (127), t = 0, we find (<^,^) = (^q— 5)o, whence (151) follows. 
The reader will notice that (151) can also be written in the form 
of a Stieltjes-integral 

<l>o=-Trfit)d<!>(t). (152) 

Jo JO 

It is, of course, easy to obtain Theorem X from (152). Applica- 
tion to (123) yields, according to ?o=/o> 

which together with (139') proves the theorem. 



46 


SOLUTION OF PROBLEMS I AND II 


Application to the boundary temperature in Problems I and II , 
In the case of Problem la we have 

2II,{0) = E,(0) = 1 

1 


yielding thus 


/o = 


V3 


.(153) 


for the solution mentioned in Theorem IV. According to (84) we 
therefore get 




.(154) 


Under local thermodynamical equilibrium we thus find, according 
to (98), the value /o 

7T 4 _ V rp 4 

-*■ 0 ^ e 

for the surface temperature of model la, independent of the 
variation of the absorption coefficient within the slab. 

Let us, now, find a formula for the boundary temperature of 
the model II a. In order to apply Theorem TX we first suppose 
that € ^ 0 everywhere. We find for the iV'-solution B (r) of (130) 

^o-"| r-f’dt = f + r^-q'dt. 

^ lO /o J 0 ^ /o J 0 Ot fo J 0 ^ 

According to (132) and (153), the first term on the right equals 

4^0- 

From Theorem IX we therefore obtain, on account of Fq = Fq + , 




I f -q'dt, (L55) 

o'* Jo a 


for the determination of the surface temperature of model II. 

It should be noticed that (155) holds also in the general case, 
6^0, provided that 




^ I pdx, 


integrated through the slab, is finite. In the next section we shall 
prove that q{T) increases with increasing r. Anticipating this 
result, we infer from (155) that, in the case e ^ 0 everywhere, the 
boundary temperature is greater than in strict radiative equili- 
brium, 6 = 0, the surface flux ttFq being prescribed. 
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This minimal property of the boundary temperature in the 
case of strict radiative equilibrium implies, conversely, the in- 
creasing of the remainder function (/(r). Otherwise the energy 
liberated could evidently be distributed in such a way that that 
property becomes violated. 


§15. PROOF THAT q (t) INCREASES 

Let us, for a moment, reconsider the hypotheses made in the 
general case treated in §§ 11-14. The hypothesis (112) concerning 
the behaviour of H (.r) for small x was only required for the proof 
of formula (114), of Lemma 4 and of Theorem XI (existence of 
/'(r)). 

For the proof of the Theorems VI-X we needed, however, only 
the suppositions stated in (116), (117) and (119), the latter being 
of particular importance because it is responsible for the asymp- 
totically linear character of the solution. 

Let us, now, throughout this section, suppose that H (a;) is of 
the form 

e~^'^dp{s), (156) 


p (s) increasing with increasing s, (119) is satisfied if and only if 


From 





1 

2 ‘ 


1-/1 


dp (s), 


(156') 


we see that ( 1 17) is certainly fulfilled, with a = 1. The last integral 
is readily seen to converge also for n^O, making thus the formula 
valid for all 0. 

On applying Schwartz’s inequality to in the same way as 
we have done in deriving (76), we find 


for all n, implying thus 


dx 'H„ {z) ■ 


(157) 
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The analogue of (78) subsists therefore in the present general 
case, the proof being the same, 


.(158) 


d II Ax -a) ^ 

-y--';/ — <0; x>a. 
dx {x) 

We now begin with the proof that, under the hypothesis (156), 
(156'), the remainder function q (r) of Theorem VI increases with 
increasing r. 

Transforming the integral operator A by partial integration, 
we find the identity 

(t) — A {u)^ = Uq II 2 (t) -h J II 2 (t — t) du {t) — J //g {t — t) du (t), 

(159) 


Stieltjes-integrals being employed because we do not wish to 
make any use of the differentiability of/(T). (123) can therefore 
be written 


Let, in the sequel, u{t) denote any positive and nowhere 
decreasing function that satisfies the integral inequality 

+ (161) 

for all T. Continuity of u is not required. An example of a function 
u{r) is furnished by 1, for in (121) we have H 2 > The set 
of values taken by all functions u (r) at a particular point r has 
a greatest lower bound U—U(t). u^O and (161) implies 
U (r) ^ /fg (r) > 0. Furthermore, (7 (r) is evidently nowhere 
decreasing. From (161) and from u^U we see that 
u(r)^A{U), + H 2 {T) 

holds for all functions u, and therefore for their greatest lower 
bound U. U (r) is thus itself a function u, in fact, the smallest 
function of that class. 


Lemma. U is continuous for r ^ 0. In the inequality 

U{t)^A{u), + H2(t) (162) 

the equality sign takes place at all points not being inner points 
of a constancy interval of C7. t = 0 is included herein. 
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We postpone the lengthy proof of this main lemma to the end of 
the section and show first that U cannot have a constancy in- 
terval, thereby proving that U increases and satisfies the same 
integral equation as (/(t). This would, however, imply U = q, 
since a bounded solution of the homogeneous equation vanishes, 
according to Lemma 2. Supposing that U have a constancy 
interval we find, according to the lemma, from (159) 

and (162) 

= a, 

> a, 

( 163 ) 

T lying between a and j3. This can, however, be shown to lead to 
contradiction. According to (157), the second term decreases 
with increasing r. According to (158), and to dU^O, the third 
term certainly never increases. Furthermore, H 2 (t — r)IIl 2 (r) 
increases with increasing t, the denominator decreasing and the 
numerator increasing. The fourth term in (163) is therefore also 
a nowhere increasing function of r. Altogether we find that the 
whole left-hand side of (163) decreases, a<T<j8, in obvious 
contradiction to the right-hand side. 

Proof of the lemma. For a nowhere deereasing function U the 
limits U (t — O) and (7 (r 4- 0) exist always. Since U is the smallest 
of all functions u, we find 

U(r^0)=U(T), 

since otherwise u(t) = U (t — 0) ^ ?7 (r) would be a smaller func- 
tion u. 

In order to prove the continuity of U we must show that 
1 J = 0, where 2ri = U {t^ + 0)-U (164) 

The right-hand side of (162) being a continuous function, we 
derive from C7 (t) ^ i7 (tq -f 0), t > tq , 

U (T)^A(i7)^-f ^3(t)4-77, To<T<To-fS, (165) 

with a suitable 8 > 0. On introducing the auxiliary function 


“ Ihir) 


lh{r) 




1, T in (to,To + 8), 

0, T not in (to,to + S), 


(166) 


H M P 


4 
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we first find the obvious inequalities 

t in {to,to + 8), 

A(A)t, t not in (tq.to + S). 
It is geometrically evident that the function 


h(r)< 


.( 167 ) 


u (t) = U (r) — 7]h (t) 

never decreases and that Furthermore, the second in- 

equality (167) shows that, outside the interval (To,To-h8), u 
satisfies (161). On the other hand, we infer from (165) and from 
the first inequality (167) that (161) is also satisfied by u within 
that interval, u is therefore a member of the class introduced 
above, and u ^ f/, which obviously implies t; = 0. 

It is immediately seen that, in (162), the equality sign must 
hold at T = 0. Otherwise we could make the value of U (0) slightly 
smaller without affecting (162), in contradiction to the definition 
of the smallest ^-function. 

A point T shall be called a ‘proper point’ of U {r) when 
U (t') > U (t) holds for every r' >t. According to the continuity 
of U, an ‘improper point’ is thus either a left end point or an 
inner point of a constancy interval of C/. A point that does not 
lie inside of such an interval is therefore either proper or a left end 
point of such an interval. In the latter case, and for r > 0, the point 
is certainly a right-hand limit point of proper points. According to 
the continuity of U it is thus sufficient to prove the equality sign in 
(162) for all proper points tq . We define a number ^ 0 by setting 

(169) 

and show that ry > 0 is impossible. As U is continuous, (165) must 
hold in a sufficiently small interval (tq, tq-I-S). This time we use 
another auxiliary function 

7 . V |t^(T)-i7(To), T in (to,To + S), 

^ 1 0, T not in (to,To + 8). 

According to 0<C7^1 we have t]<1, A(t) is again seen to 
obey the inequalities (167). The function u, defined by (168), 
belongs then again to the above class, since u satisfies (161) and 
decreases nowhere (geometrical evidence). This completes the 
proof of the lemma. 
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§16. OTHER PROPERTIES OF q (r) 

The equation (96) can be immediately generalized, 

This equation can be checked by direct differentiation, taking 
account of (123). There should be an additional integration con- 
stant on the right. On comparing, for r = 0, with (139) we see, 
however, that it has the value zero. The left-hand side of (170) 
can easily be transformed by partial integration, thus yielding 


/: 


H^(\r-t\)dq(t) = II ^ (r) - qo (r) (170') 


From (157) and from follows the existence of the limit 

A.lim al, 

r=oo /H-1 V*^/ 

thus yielding (a:) = where S = S(a;)->0 as x->co. 

According to (156) we have 

Hn+l (^) > J (s) > I 8~^^dp (8), 

p increasing and € being arbitrary, whence follows A^l. This 
implies X=l 

From (158) we infer, furthermore, the existence of 

H„{x-y) 


lim 




x=« //„(^) 

From we obtain, by differentiation of this 

limit relation with respect to y, (y) = <j>(y). This implies, together 
with <^(0) = 1, <^(y) = c^. According to (157) and (158), we find 
at the same time that 

Hn (t - <) > (t). 

Since, now, q{t) increases, we find from (170') 




thus yielding 


J e^dq{t)‘^l — 




.(171) 


4-2 
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Let US, in (156), make an additional hypothesis concerning the 
behaviour of p (s) at s = 1. We suppose that 

' dp {s) 

II 8—1 




- =CX). 


.(172) 


When p has a continuous derivative this is certainly true if 
/)'(0)>0. This takes place, for instance, in Milne’s case, 
p (s) = -Hog 5. Under this hypothesis we always have 


i: 


e‘dq{t) = l-qo. 


.(173) 


It is remarkable that this general relation does not contain H or 
p explicitly. 

Proof of (173), We first show that the functions 


have the properties 




(«) ^‘’ds 


.(174) 


.(174') 


According to 1 as 5->co, the second property is seen 

to be a consequence of the first one. Interchanging the order of 
integration wo find 


(*oo rcoffi—u 

HJx)e-dx=j^ —dp{s), 


which is infinite according to (172). 
We need the inequality (c = (0)) 


J. 


C-^3(|T-«|)dr^ 


e^{li{x)-\rc}, t^x, 


ce^, 


t>x. 


.(174") 


For t > X the left-hand side equals, indeed, 

I H^(t — T)dT <e^ | H^(t — T)dT<ce^<ceK 

Jo Jo 

For t^x, however, it equals 

(^ — t)(ZtH- H^{T — t)dr, 
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The first term is smaller than while the second term becomes, 

on substituting t = ^ + 5, 


A ^ 

eM e^H^{s)d8, 

being obviously less than (.r), which proves (174"). According 
to (174") we get, on multiplying (170') by and integrating 
through 0<T<x, 




f*.r /•« 

} ^dq + 

Jo Jr 


e^dq. 


Dividing through by l^{x) and proceeding to the limit x->co, 
leads, according to (174'), to 



whence, with respect to (171), (173) follows. 

By partial integration, (173) can be transformed into 

(173') 

Numerical remarks on Milne's case. Let us first find the value 
of the constant a in the formula (97). The limit of the left-hand 
side as T->oo evidently equals the limit of 

f q{t)E^{\T-t\)dt, 

J (ii 

(X) being an arbitrarily fixed quantity. This integral is, now, 
included between the limits 


q{ci)) r J?3 (I T — ^ |)d^, q^ f 
J cj J a) 

Here, the integral tends to | as t->oo, whence 


and according to the arbitrariness of w, 

a=5'7». (175) 

From (97), t = 0, and from (175), 

{q,E^)=lq^-h ( 176 ) 


Furthermore, from (96), t = 0, or from (139), H = 

{q>E^=\. 


( 177 ) 
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Practical computation of a number means including it between 
sufficiently narrow limits. This can sometimes be done by means 
of mathematical artifices. The two relations (176), (177), to- 
gether with the fact that q increases, suffice to find the value of 
up to 1 per cent., 1 1 (178) 

or, in rounding off the second digit, = 

We have, in fact, from (176) and (177), 

which, according to 2E^ > for r > 0, implies 

?oo < f?oo ( 1 > ^^ 3 “ ^^^2) + 8 = + h 

i.e. the right-hand inequality of (178). On the other hand, from 
3E^>2E^, ^00 

(-^4”" 3-^3)^? ^ 

Since here the integrand vanishes at ^ = 0 and at ^ = 00 , we obtain 
by partial integration ^ 

Together with (176) and (177), this implies the second inequality 
of (178). 

We now consider the law of darkening (101) and notice that the 
second term in the parenthesis decreases with increasing 6, 
0<7r/2. Milne’s first approximation consists in replacing J by a 
constant such that the emergent radiation gives the required net 
flux ttF, thus yielding the value f for the constant. F being a 
definite average of the emergent radiation, we see that the 
second term lies, for 9 near 7r/2, below Milne’s approximation f 
and rises, as 9 decreases, above it. The smallest value of the second 
term equals = l/\/3, while the greatest value (q, surpasses f . 
We may add that | < ^ 

Wo have, in fact, from (177) and from E^ < e~*, 

{q, e-') = {q, e-f-E 2 ) + i<^^ + l< 0-69. 
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§17. PROBLEM Tfo FOR PARALLEL 
INCIDENT RADIATION 

Parallel radiation carrying finite energy is to be considered as a 
limiting case. Let 

r' = {7r^e\n 

be the direction of an incident ray, and let Aa> be a solid angle of 
directions containing r\ We set 

[0, r not 111 Aa>, 

for the incident radiation. When Aa> shrinks down to the single 
direction r\ we speak of parallel radiation of the direction r' and 
of the flux 

through the unit area normal to r'. On proceeding to the limit 
Aa>->0, Milne’s fundamental integral equation (53) becomes 


J5 (t) = A (B)^ + 5 = sec 0'. 


.(180) 


We define (180) as the integral equation for parallel incident 
radiation of direction r'. 

According to Theorem VII, (180) has a positive solution. Let 
us, in particular, consider the iV-solution, i.e. the smallest 
positive solution of (180). We put 

B(T) = Sg,(r), (181) 

gg being the V-solution of (180) for S=l, 


Theorem XII. The A-solution of (180) can be explicitly ex- 
pressed in terms of the solution /(t) = r -f g (r) of Problem la. 


9s (t) = f if, e-«0 1 ~ - f (<) dt 
\ ^ Jo 

Proof. (182) can also be written in the form 
9s {r) = {/. e-0 { {t) d<}. 


.(182) 


( 182 ') 
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First we show that the function 

(182") 

solves (180) with a suitable S. From (182"), 

sh + h^=f, hii)) = 0. (183) 


Formula (114) is applicable, h being continuously differentiable 
and satisfying In! = 0 (r) for t large. On setting 

p(t) = A(t)-A(A)^, 
we find from (114) and from (183), 

•5p + p'=/-A(/) = 0. 

Hence p (r) = const. Now we show that in 

= + ( 184 ) 

a equals J. It is obvious from (182') that (r) is a bounded func- 
tion of T. The boundedness clearly implies that it is the A’-solution 
of (184). Applying, now, formula (152) to (184), 0 = ^6“^^ we find 


Jo 

If (7, (0) is computed from (182), we obtain, according to/o = l/\/8, 
a = }. This completes the proof of our theorem. 

We note the special cases t = 0 , oo, in (182), 


<78(O)="^s(/,e-*0, 




thus yielding 


9s W 

<7,(00) V3’ 


(185) 

(186) 


(185) being expressed in terms of the solution /of the Problem la. 
In particular we have the formula 


for the ratio of the limiting temperatures. 

The fundamental limit relation (56), J = B, Foo = F, is easily 
seen to remain vahd in the case of parallel incident radiation. 
According to the boundedness of B{t) the right-hand integral 



SOLUTION OF PROBLEMS I AND II 57 

represents also a bounded function of r, whence = 0. We may 
collect the above results in the 

Theorem XIII. The iV'-solution B(t) = Sg^ir) of (180) corre- 
sponds to vanishing net flux, F = 0 , i.e. it gives the solution of the 
Problem 16 (insolation problem) for parallel incident radiation of 
direction r' = (77 — 0 \ </)'), s = sec 0 \ 

By partial integration we find from (182'), according to (185), 
that (T)-g^ ( 00 ) equals 

A(T)=f e»'dg'(<)-(l-(Zo) + l-7. 

Jo ^ 

up to a positive factor. For we have, according to (173), 
A < 0 for all r, whereas for s> I, 

A(aD)>l-^>0. 

Y o s s 

A(t) being an increasing function, we infer that the equation 
A(t) = 0 has at most one root. Furthermore, on differentiating 
(182'), and using (173), we find that, for 5^1, g,^ (r) increases with 
increasing r. Altogether we note the following behaviour of the 
solution B(t) of (180): 

B(r) increases. 

\/3 >s> 1: B(t) lies first below the value B^^ and then rises 
above it. 

s > \/3: B (r) lies entirely above B^ . 

5 ^ 1 (5 = 1 corresponds to normal incidence) is the only case 
where B is monotonic. 

Oeneralization, Model 16 with parallel incident radiation, and 
with F = 0, mainly applies to the upper layers of a planetary 
atmosphere. The hotter solar radiation will actually have a 
smaller absorption coefficient than the cooler radiation of the 
atmosphere. The next idealizing step would thus consist in 
making the former coefficient a constant fraction n of the latter 
coefficient a. 
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In order to find the integral equation we must go back to the 
fundamental flux equation (54). The first term on the right is the 
only one that contains the radiation from outside. In this term, 
therefore, r must be replaced by ut. Proceeding as before to the 
limiting case of parallel radiation, we find, J—B, 


:aScos 0 'c-^^«®^^' + 2 BE^(i-T)dt -2 BE^(T-t)dt, 


= const. On differentiating this relation and on dividing 
through by 4, we find 

7} Si 

B(r) = A (B)^ + ~ 


Since this is of the same type as (180), all formulae of this section 
remain unchanged, with the only difference that S is to be 
replaced byn/Sand that 5 is now sec O'. We have, instead of (181) 

and (182), ^ sec0' 

B{T) = nSg,,,Mr), £ = 


§ 18 . THE EMERGENT LIGHT 
It will be convenient now to introduce the new notation 


r,(a) = (7 fV-^(/,(0d^ a>(c7) = or re-^f(t)dt (188) 

Jo Jo 

The emergent radiation (law of darkening) becomes then 

Problem la: / (0, 0) = (sec 0), (189) 

and in the case of parallel incidence. 

Problem 16: I (0,6) = ST 6), (189') 


From (183) we find, by partial integration, 

I e~^^h'dt = (7 1 er^^hdt, 
Jo Jo 

thus yielding, again with respect to (183), 

^00 

O ((j) = <7 (5 + (7) J e~°^hdt. 
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This gives, together with (182) and (182”), the formula 



■i 5 + O’ 

We have therefore the 


Theorem XIV. In Problem I ft for parallel radiation with the 
angle of incidence O', the emergent radiation can be explicitly 
expressed in terms of the same radiation for Problem la, 


i(o,e) = is 


cos fl'O (sec 9) O (sec O') 
cos 0 + cos 


(191) 


In the general Problem I, with ]mrallel incident radiation, but 
with an arbitrary flux ttF ^ 0, the emergent radiation is, of course. 


1(0,9) 


_3 "■ 
“4 __ 


F + S 


cos 9'(S> (sec 9') 
cos 9 H- cos 


?«'n® 


(sec 0). 


This applies to binary stars, where the incident radiation is due 
to the other comj)onent. 


§19. PROBLEM 16 FOR ARBITRARY 
INCIDENT RADIATION 

Since arbitrary incident radiation can be obtained by super- 
position of parallel bundles of different directions, we should be 
able to express the solution of Problem I ft, for any given incident 
radiation, in terms of the solution for parallel incidence. For 
given incident radiation I (0, /) > 0, a solution of (53), € = 0, J = 5, 
is, in fact, obviously given by 

IT 

1 /• 27r /• 2 

= - /(O,7r-0',^')f7scc0' Wsin0'd0'cZ<^' (192) 

^Jo Jo 

Let us show that this solution corresponds to vanishing net flux, 
F = 0. We note that, for any value of r, (r) decreases with in- 
creasing s. This is physically evident, for it means that the 
temperature becomes smaller when the incident radiation is more 
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oblique. Remember, in fact, that g^fj) is the iV'-solution of (180), 
5=1. The above assertion follows simply from the fact that the 
second term on the right of ( 1 80) decreases with increasing 5,5^1. 
We have, in particular, g^ (r) g g^ (r), and from (192) 

l(0,r')dco'. 

The boundedness of -B (t) implies, as before, the vanishing of the 
net flux. 



CHAPTER III 

DISCUSSION OP PROBLEMS III AND IV 

§20. SOLUTION OF PKOBLFM III 

Problem III has been analytically set up in §7. We remember 
that the flux ttF is constant and that the limit relation (44) holds, 
in consequence of the hypotheses, that the law of scattering is 
symmetric in the directions of the incident and scattered ray 
(11'), and that the same amount of radiation is scattered in 
opposite directions (IT'). 

Theorem XV. Problem III always has a solution 

(193) 

with the property 

f(r,r) = T + q(r,r); 0<q<l, (194) 

the limits thus being independent of the law of scattering. 
Under the hypothesis 

y(T;r, r') > const. >0 
the solution is unique. 

Proof, The reasoning is essentially the same as in Problem la, 
i.e. in the special case of uniform scattering, y = 1. We put 

= e7(T,r) = T+l, (195) 

and denote by the respective intensities found from (47) 

and (57) in setting J — J = J, On inserting these intensities 
into the right-hand side of (41) we obtain the expressions 
A(Ji),,„ A(J\„ 

respectively. 

Using the more convenient relations (99), (99'), instead of (47) 
and (57), we readily find 

J , (T-I-COS0, d<7r/2y 

1^ [T + C0sfl+ I COS0 I 0>7t/2, 

T-f-l+cos^, 0<7r/2, 

T-f 1 + COS0 — (1 -f cos0)e"“'^>^®‘=^i, 0>7r/2. 

These equations imply, in particular, the important inequalities 


and / (r, r) = 


itions imply, in particular, the important i 
(t, r) ^ t + cos 0, 7 (t, r) ^ T -h 1 + cos 0, 
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for all directions r = (0, <f>). Inserting this into (41) and taking into 
account the relation 

J cos 0'y (t; /, r) dco' = 0, 
we obtain the important inequalities 

(.'T, r). 

According to the positivity of the operator A, we can, now, apply 
the same reasoning as in § 9. On setting 

{*^n)T»ri 

we infer that lies above J ,^ , but that all lie below J. The 
limit function / (r, r) is seen to satisfy the integral equation 

/(T,r) = A (/),,„ (196) 

and to lie between the indicated limits t and t+ 1. 

We set J (r, r) = c/(t, r) 

for the solution of J of Problem III and determine the constant 
c as in § 9. The intensity obviously lies between the limits 

c (t + cos0) < I { r , r )< c { T -{- 1 + cos0), (1^7) 


for all directions r, whence 


hm = c 


uniformly for all directions r. This gives, according to the 
fundamental relation (44), 




(198) 


The proof of the uniqueness of the solution follows the same 
lines as in § 10. It is sufficient to repeat the main part of the proof, 
i.e. to show that a solution J of (196) with the properties 


T,->00, 

V=ao J 


Ty being a given sequence of numbers tending to infinity, vanishes 
identically. 

According to the hypothesis made about y, and according to 
(41), we have 

J = ^ J/ydoi' ^ const. J/ cos^ O'do)'. 
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On comparing this with (44) and (199) we get i^ = i^ao = 9. The 
incident radiation (at t = 0) being zero we thus have, for r = 0, 

F = = - { / (0, r) cos 6dct} = 0y 

which is possible only if the emergent radiation vanishes. Thence, 
and from (99), t = 0, wo infer that the Ergiebigkeit must vanish 
too, q.e.d. 


Remarks on the uniqueness. When the hypothesis made about 
the law of scattering is abandoned. Problem III could have 
several solutions. This takes place, for instance, in the easily 
integrable limit case, where half the radiation is scattered back- 
wards while the other half goes in the original direction. Under 
the hypothesis, however, that 


lim 

T=00 


I ('T, r) 


holds uniformly for all directions, i.e. that the radiation becomes 
isotropic at great depth, the solution is readily found to be unique. 
For we find from (41), (44) and (194) that then 


r=00 y(T,r) 


holds uniformly for all directions. As in the proof of Lemma 2, 
§ 12, we infer, from the positivity of the operator, that J = JF/* 


§21. DISCUSSION OF PROBLEM IV 
IN THE CASE ((35) 

The optical thickness r of the slab being finite, we consider the 
case (65) of constant radiation incident at the inner face t = r^. 

For the radiation coming from above we use again (99'), while 
the radiation from below is given by (61), or by the more con- 
venient formula 

e{T* — T)SGGd 

+J e~^ J {r-{-scosd,r)dsi 0<7r/2 (200) 
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Theorem XVI. Problem IV, with the boundary condition 
(65), has a unique solution. The Ergiebigkeit lies between the 


limits 




.( 201 ) 


the limits thus being independent of the law of scattering. 


Proof, This theorem has, in the case of uniform scattering, 
been proved by Schwarzschild. 

The integral equation (63) of Problem IV can be written in 


the form 

J{T,r) = L{J),^, + H(r,r), 


where 

H(T,r) = ~ f 

....( 202 ) 

We set 


....(203) 


In order to compute the quantities 

L (Ji)r,r+ K (r, r), L + K (r, r), 

we first find, from (91)') and ( 200 ), the intensities and I corre- 
sponding to the Ergiebigkeiten and «/, and then insert them 
into the right-hand side of (41). After some simple computations 
we get 

/* fr -f cos 0 -{- ( 1 — cos 6 ) e-(^*-r)sec 9^ 0 ^ 

T-f-cos0-f 1 COS0 I 0>7r/2, 




and 


I (T,r) = 


^+i' 


T + 1 + COS 6 — cos 0 e-<T*-T) 8 ec 0 ^ 0 ^^ 2 , 

T+ 1 + COS0 — (l + coS0)e~’''®®®®', 0>rr/2. 


These relations imply, for all directions r, the important in- 


equalities 




■ + COS 9 
T*+l ’ 




On inserting the intensities into the right-hand integral of (41), 
we obtain, according to 

jeos 0 'y (r; r', r) da>' = 0 , 
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the final inequalities 

L ( + H (r, r) > I* (r, r), (204) 

L (/),., + II (r, r) < I* = J (r, r) (205) 

According to the positivity of the operator L, we see as before 
that there is a solution between the indicated limits. On setting 

(r, r) = L {J,X.r + U (t, r) (206) 

Ji being defined by (203), we can write (204) in the form J^> 
This implies L{J^>L{J^) and, according to (206), 7 i=l, 2 , 
J 3 > Jgj generally J Furthermore, the relation 

J — J,i+i — ^ n) 

and the inequality /— > 0 show that generally < J- The 

limit function J of the functions thus exists and has all the 
properties indicated in Theorem XVI. J is, moreover, seen to 
be the iV'-solution of (63). 

In the present case of a finite r^, the uniqueness is inferred by 
a classical conclusion. First we note that 


’ J _ *1^ J + 


in particular, that L( 1 )^,.^A <1 


holds with 'a suitable constant A. The difference D (r, r) of two 
solutions of (63) satisfies the homogeneous equation D = L{D), 
Denoting by M the maximum value of | D | , we find 


\D\^L{\D\)^ML{l)^\M\ 


in particular M ^ XM, which, according to A < 1 , gives Jf = 0 . 


Large optical thickness. In the case of uniform scattering, 
Schwarzschild recognized, that for t^->oo, model IV goes over 
into model III. This holds, more generally, for any given law of 
scattering, of the types ( 11 ') and ( 11 "). 


H M P 


5 
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Theorem XVII. For = c being fixed, the Ergiebigkeit 
J (r, r) tends towards cf (r, r) as the optical thickness r* increases 
indefinitely. 

The main part of the proof is to show that the Ergiebigkeit 
J (r, r) tends towards a limit function as -> oo. We may content 
ourselves with a brief sketch of the proof. From the integral 
equation (63) it can be first recognized that the Ergiebigkeit is 
uniformly continuous for all r*. Arzela’s selection theorem then 
guarantees the existence of one or several limit functions. Such 
a limit function necessarily lies between ct and c (r + 1 ) . From (11) 
and (202) we have 

H (t, r) < I* = CT* c^ -> 0, 

as T*->oo. On proceeding, in (63), to the limit t* = oo, along any 
special sequence of numbers r*, we find that every limit function 
must satisfy the integral equation of Problem III. This equation 
having, however, only one solution cr for r large, the Ergiebigkeit 
can have only one limit function, as r*->oo, namely c/(t, r). 

§22. NARROW LIMITS FOR THE SOLUTION 
IN SCHVVARZSCHILD’S CASE, 

In the case of uniform scattering the integral equation becomes 
(64), (64'), the second term on the right being (66). On setting 


J(r) = l*u{r), 
u becomes the solution of 

u{r)==L {u)^ + \E^ (r* - r) (208) 

Let V (t) denote the solution of 

= (209) 

u is then the Ergiebigkeit in the case 

/*=!, (208') 

while V corresponds to the case 

/ (r*, 0) = cos 0 ; e< 7r/2. (209') 


On introducing, for any function, the general notation 

O (t) = $ (t* — r), 
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we readily find that the operator L, given by (64), satisfies 

L (0\ = L mr*~r= LWfr (210) 

From (207), y = 1, we have, according to (49), w = 2, 

l = i(l)^+P2+P2. (211) 

Furthermore, we find by means of (79), (81), (82), n=\, 

t=L (t), + + 4 (^3 - Fj) (212) 

(208), (210), (211) imply 

u + ll—l = L(u + ll—l), 

whence, according to the uniqueness property proved in the 
preceding section, m + I?=1. (213) 


(213) is physically obvious. For u is the Ergiebigkeit when the 
radiation incident is zero at r = r* and one at r = 0. -f- w is there- 

fore the Ergiebigkeit when the radiation incident is one at both 
faces, which yields (213). 

From (209), (211), (212), 

v — — t — (214) 

In the sequel we make use of the 

Lemma. (t^, d) < /g (r*, ^), 6 < 7r/2, implies (r) < (t) for 

the corresponding Ergiebigkeiten. 

This follows immediately from the fact that the Ergiebigkeit 
is the iV-solution of (64'). 

Narrow limits for u (r) can be found by comparison with the 
solution of Problem la. On setting 

qi = q{r*), q^^ j ^e-^q{T* + a)d3, (215) 

we infer that and that qi and increase with increasing r^, 

having both the limit q^o as r* -> oo. In the model la, the intensity 
of the radiation coming from below is given by (100). Let us in 
the case F = ^ denote this intensity by I, Since q (r) increases, we 
then find from (100) 

T* + + cos 0<I (r*, 9)<t* + q^ 4- cos 0, (216) 

for 0 < 7r/2. Taking I as the radiation incident at the inner face 


5-2 



68 


DISCUSSION OF PROBLEMS III AND IV 


T = T* of model la, we first see that the corresponding Ergiebig- 
keit is precisely f{T) = T + q{T) (according to the theory of model 
la). On the other hand, the Ergiebigkeit turns out to be 
{T* + qi)U'\-v resp. (t* -\- q^u-{- v, 
when the radiation incident at t = t* is given by the left respective 
right-hand term in (216). The lemma shows therefore that, for 
all T between 0 and r*, 

{r* + qi)u-\-v<f<{T'^-\-q^)u-\-V, 
which, for the function 

f-f=2T-r* + q-q, 
implies the inequalities 

/-/> (T* + qi)u + v-{r* + q2)u-v 
and /— /< (r"^ + ^ 2 ) + qj) u — v. 

When, here, u and v — v are expressed in terms of u by means of 
(213) and (214), u is found to lie between the limits 


T* + <?1 + </2 

■’■* + ?1 + 

(217) 

By means of the notation 



It 


....(215') 


we can give (217) the form 

J I I . r* 

^ - V* + 2m(T*) I T* + 2m(T-*)‘ 

(217') 

The use of (217) requires, of course, a numerical knowledge of the 
function qir) of Problem la. The limits given in (217) are much 
narrower than SchwarzschikVs limits indicated in Theorem XVI, 
since q(T) is known to lie between the narrow limits 0*577 and 
0*71. The factor ^(t’*) on the right of (217') is about 0*05 for 
T^ = 0 and decreases very rapidly as t^->oo. 

If, in (217), q is replaced by Milne’s approximate value f, we 
obtain Milne’s approximate form of J, 

T*l±l 

T*+r 
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while for large, the true approximate expression is 

r^, r + q{T) 

In the case of uniform scattering, Theorem XVII is an evident 
consequence of (217) or (217'). 

§23. THE NET FLUX IN 
SCHVVARZSCHILD’S PROBLEM 

The flux constant F is not given in Schwarzschild’s problem. For 
astrophysical purposes it is, however, of importance to know 
as a function of r'^. We shall now include this function 
between narrow limits. 

First we must find the analogue of the equation (54) in 
Schwarzschild’s model. The first term on the right must, of course, 
be omitted, since the radiation incident at the outer face is zero. 
Instead of this term there is the analogous term for the radiation 
incident at the inner face, 

Q(r) = - f cos0c-(^*~'^)^^^^/(r*,r)da> (218) 

J + 

Furthermore, the upper limit of integration in (54) is to be 
replaced by r^. We thus obtain Milne’s flux equation 

IF = J (t) {t -T)dt-^^ J (t) E^(T-t)dt+lQ ir) ■ 

(219) 

For r = 0, T*, we have 

lF = j^^JE^dt + ^Q{0), ( 220 ) 

IF=- r jE^dt +IQ (t*), (220') 

* Jo 

whence f {J -J) (220") 

Jo 

Let P resp. F' denote the flux constant in the case (208') resp. 
(209'). Taking account of (49), 7i = 3, 4, we note the following 


combinations in those two cases: 

Q = 2F3(t*-t), F = P (221) 

and J = t;, Q — 2E^{r* — T), F = F' (222) 
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In the present case (65) we have 

F = I*P 

for the net flux. 

Let us, first, find the relation between and P. On applying 
(220") to the case (222) we find, according to (214), 

f tE^dt — T* r uE ^dt \ + E . 

Jo Jo 

The first integral is readily computed from (80), n = 2, while the 
second one can, by means of (220), taken in the case of (221), be 
expressed in terms of P. Altogether, the simple relation 

= (223) 

is obtained. ^ 


The following remark will be of use. The inequality 
A (t*, 6 ) < 1 2 (t* 6 ), 6 < 77/2, 

implies the inequality F^ < F^ for corresponding flux constants. 
This follows from the lemma of the preceding section and from 
(218), (220). 

We now recall the inequality (216), where the incident radia- 
tion I gives precisely the net flux F=^^. Taking account of the 
above remark and of (221), (222), we find therefore 

{T*-\-qj)P + F' <^<{r* + q2)F+F'. 


If this is combined with (223), we obtain simple limits for 




4 7* 4 7^ 


(224) 


By means of (215') we can write this 



This approximate expression of IjP is, even for small t^, very 
accurate. Its value for t* = 0 turns out to be 0-95, differing only 
5 per cent, from the true value IjP—l, 

In comparison hereto, we note that the right-hand sideof (224') 
gives 8 per cent, at t* = 0. For larger t* the accuracy is extra- 
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ordinarily good. According to m->q^ as we find, in 

particular, from (224'), 

(225) 

F 


If q is replaced by its Milne approximation we obtain Milne’s 
approximate formula 


The accuracy is inferred from the fact that 3m(T*)/2 varies 
between the narrow limits 0*95 and 1*06. 



CHAPTER IV 


EXPLICIT SOLUTION OF CERTAIN 
INTEGRAL EQUATIONS 


§24. THE CHARACTERISTIC EQUATION 


We shall now generally treat the integral equation 

f{x)=^^H{\x-y\)f(y)dy (220) 

by means of Fourier integrals. The solution will be obtained in 
the form of explicit integral formulae. 

Concerning the kernel, we suppose that, for a certain « > 0, 
is quadratically integrable within 0<a;<oo. Without 
limitation of the generality we may suppose that 


j: 


[H {x)e^’^]^dx<QO\ s<l. 


.(227) 


This hypothesis is certainly fulfilled in Milne’s case. Quadratic 
integrability is introduced in order to apply the Plancherel 
Theory of Fourier Integrals. 

The solutions of the analogous equation 


f{x)=( II{\x-y\)f{y)dy (226') 

J —00 

are of much simpler form, being aggregates of exponential func- 
tions. If u* denote an w-fold root of the ‘ characteristic equation’ 

f +00 

l = /c(w)=j H {\x\)e^^dt, (228) 

J —00 

the function Qn-i (^) 

is easily seen to represent a solution of (226'), being an 

arbitrary polynomial of degree not greater than n—l. We must, 
of course, be sure that all integrals involved converge. This is, 
according to (227), obviously true when the real part of u lies 
between the limits - i < (^4*) < 1 

We should expect that the characteristic equation plays also an 
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important role in the theory of (226), and that its solutions show 
the same behaviour for large x as certain solutions of (226'). 

We first compute the characteristic function k (u) in the im- 
portant case, where H is of the form (156), with p (s) increasing. 
This includes Milne’s case with p = ilogs. k(u) is readily seen 
to equal 1*00 n 1 1 “I 

Lf 

In order that (227) be fulfilled, we need only suppose that 
K (0) = 2^2 finite. In Milne’s case we have 

= (230) 


the logarithm being zero for u = 0. 

Under the hypothesis (227), the complex function k(u) is 


certainly holomorphic in the strip 

|iJ('i4)|<l. (231) 

Two obvious properties are 

K{u) = f<(u), k(u) = k( — u), (232) 


showing that k (u) is real on the imaginary as well as on the real 
axis. 

(227) implies absolute integrability 


i: 


I H (x)e^'^ ldx<oo; s< 1, 


.(233) 


as seen on applying Schwarz’s inequality to (233), the integrand 
being / \ 

\H(x)e^ Je^ . 

By simple substitution x = x' + Trjt, k(u) = k( 8’\- it) can be written 
in the form 

2k{u)= j {\ — ^ e'^^^dx^ 


whence 
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If H is continuous, this inequality shows that 

/c(5 + i^)->0 as |^|->oo (234) 

holds uniformly in every partial strip | 5 | 1. Theorems of 

Lebesgue, however, show that this still holds under mere measur- 
ability of H, As a consequence we note that the characteristic 
equation can have only a finite number of roots within every 
partial interval | 5 | ^ a < 1. 

More can only be said about the characteristic roots when H 
satisfies special hypotheses. In the important case, for instance, 
where H is positive, k{u) is found to increase from k(0) to 
K ( 1)^00 when u goes from 0 to 1 along the real axis (or from 0 to 
— 1). On the imaginary axis k{u) is still real and less than /c(0) 
in absolute value. The characteristic equation has therefore the 
double root u = 0 when /c (0) = 1 (235) 

two real roots of opposite sign, however, when 

/c(0)<1<k:(1), (235') 

and at least two imaginary roots when 

k(0)>1. (235") 

Besides these it could very well have other complex roots in the 
strip I 5 I < 1 . 

Let us, now, suppose H to be of the form (156) with increasing p. 
In this case, there are, within | Q(i^) | < 1, no other roots than 
u = 0 in the case of (235), than two opposite real roots in case of 
(235') and than two oi^posite imaginary roots in the case of (235"). 
In Milne’s case, for instance, all characteristic roots are thus 
exhausted by the double root u = 

which, according to (232), shows that k can only be real on the 
real or imaginary axis. Furthermore, the equation 

' ’ j. ({'+<*)■ 

proves that k decreases when u goes from 0 to -f ico or to — ico 
along the imaginary axis. The statement about the roots of /c= 1 
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is an immediate consequence of these facts. It may be added that, 
under the condition (172), i.e. /<(l) = oo, the right-hand inequality 
of (135') is automatically satisfied. 


§25. THE INTEGRAL FORMULAE FOR 
THE SOLUTION OF (226) 


We shall now construct all solutions of (226) that satisfy the 


condition 


f(x) = 0(e-'); a<l, 


.(236) 


for X large, a being an arbitrarily fixed number less than one. 


Let 




^2n 


denote the complete set of characteristic roots within the strip 
I 5 I ^ a. If necessary we enlarge a a little such that no roots lie on 
the boundaries | 5 | = a. We set 


P (u) = (u - Ui) {u ), 

P (u) being even and real on the real axis, and 

-r (u) = -- p ^ - ^y - { 1 - («)}> 


(237) 

,(238) 


T being holomorphic and free of zeros in the strip 1 5 | ^ a. More- 
over, T (u) is even and real for real u. We want r (u) in the form 
T {u) = T^{u)jT_{u), in such a way that t^\u) is holomorphic in 
the half-plane and that T_^(s + it)->l as 5 -> + 00 , while 

T-(s-\-it) is holomorphic for s^ol and satisfies T_(5-fi^)->l as 
5 -> — 00 . This obviously unique representation is, according to 
Cauchy’s formula, explicitly furnished by 


(u) = exp 


T_ (ti) = exp 




L 277* J 

-P-ico V-U J 

r--i 

C+P+i-logr(v) 1 

L 2m- J 

U/3-i» V-U J 


(230) 

...(239') 


where > a is chosen such that the strip | 5 | ^ jS contains no new 
roots. The logarithm logT(5-}-i0 is the one that vanishes as 
<-> + 00 . Observing that r(u) is never zero, and real on the 
imaginary axis, and that r (5 ± too) = 1, we see that r (w) is neces- 
sarily positive along that axis, logr (5 + it) vanishes therefore too. 
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when — We shall later see that the integrals (239) and 
(239') converge absolutely. 

On setting 

= a_{u) = {u-l)^T^{u), 

(240) 

we state the 


Theorem XVIII. The integral equation (226) has precisely n 
linearly independent solutions satisfying (236). They can be 
expressed by means of the formulae 




.(241) 


Qn-i being an arbitrary polynomial of degree less than n, and 

1 r*+ a + s 

fi^) = -TT---. (“) du; s = ^ . 

y^TTlJs—ir) ^ 


.(242) 


Corollary. The solutions have, for large x, the form 
f(x) = SQ* (x) + 0 

u* being a characteristic root within | 6' | < a, and being a 
polynomial of degree less than the multiplicity of u*. 

It must be emphasized that, while (226') has precisely 2n 
solutions with the property (236), (226) has only half the number 
of solutions. 

The abscissa of integration s in (242) can be moved without 
change of the integral, provided that no pole of (i^) = ^ (« 4- it), 
i.e. no characteristic root is met by it. / (ir) vanishes for x<0. 
Application to Milne's case. Here we have, according to (30), 




r_(u) is given by (239'), while 


U/ 


.(244) 


c being an arbitrary constant of proportionality. If (242) is 
considered as a Fourier integral, we find by the Fourier theorem 


^ (5 + 27) = f / {x) 
V 277 Jo 
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This represents <f} (u) for all 5 < 0, since the integral converges and 
is holomorphic in the left half-plane. We can obtain herefrom an 
explicit formula for the law of darkening in Problem la. From 
the last formula, ^ = 0, and from (188) and (189) we get 

I (0, 0) = IF V 27r sec 0<t> ( - sec 0) (245) 

The proportionality constant c in (243) remains to be determined 
such that the corresponding solution/(.r) is for x large. From 

(101), I ^0, = ^ which gives, together with (245), 

-y ~ = — lim scf) (,9). 

On the other hand, from (239') and from (244), we find that this 
limit equals c. We therefore obtain the law of darkening in the 
explicit form 

/ (0, 6) = F ( 1 -f cos 0) ( — sec d), (240) 

being determined by means of (239') and (243). The integral 
in (239') could be found by numerical integration. 

PlanchereVs theorem on Fourier integrals. It is convenient to 
state here the main theorems of the modern theory of Fourier 
integrals as far as we need them in the next sections. A complex 
valued function a (x) of the real variable x is called quadratically 
summable (q.s.) over —oo<x< 4-ooif 

r+oo 

j \a{x)\^dx 

J —CO 

is finite. According to Plancherel, the Fourier transform 

1 

Ait)=^~~ e^^^a{x)dx 
V27rJ -00 

exists then in the sense of mean convergence. A (t) is again q.s. 
over —oo<t< -f oo, and a (x) is conversely the Fourier transform 

oiA(x), 1 ^+co 

a lx) = A (t) dt. 

V27rJ -CO 
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Two q.s. functions a, b and their Fourier transforms A, B 
satisfy the Parseval relation 


a{x)b{x)dx==\ A 
J —00 J —00 


{x)B(x)dx. 


In particular, we have 




^dx. 


The integral (242) is to be understood in the same sense of mean 
convergence. 


§2(). PRELIMINAKIKS 

Properties of the function t { u ). We know already that t ( u ) is 
holoniorphic and different from zero within |5|^/3 (j3>a). 
Moreover, we know that r(u} converges to one as u goes to infinite 
within that strip, t ( u ) being real on the real and imaginary axes, 
we infer its positivity on these axes. The logarithm of t ( u ) con- 
sidered in ( 1 39) and ( 1 39') is therefore real on both axes and tends 
to zero as u goes to infinity within the strip, logr (u) is, further- 
more, an even function of u. 

According to (228), k {s-{-it)/V27T is the Fourier transform of 
the q.s. function e^'^H(\x\), and represents, according to Plan- 
cherel’s theorem, a q.s. function of < 1. (237) and (238) 
show therefore that logr (5 4 - it) is a q.s. function of t too. We can, 
from this fact, infer the absolute convergence of the integrals 
(239) and (239'). By means of Schwarz’s inequality we find 

47r2 1 log T_ («) 1 2 ^ I log T (w) 1 2 1 eft- 1 2 • 

J jS-icxj J P^ico I ^ 

This shows, moreover, that logT_ (u), being holomorphic ins < ^, 
is bounded in every partial half-plane 5 ^ j8' < j8. Similar facts 
are true for log (u) within the half-plane s> —p. 

The functions or+ and introduced by (240) have therefore 
the properties 

7n\u\-^^ <\a^{u) \ < M \ u\-^^ \ s^ ...(247) 
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m 1 « 1" < |<T_(w) 1 < J/| it I"; (247') 

with suitable positive constants m, J/. 

We mention, finally, that the equation t = t+/t_ can, according 
to (238) and (240), be written 

»+(w) 


\—k(u)=^P{u) 


a^(uy 


.(248) 


§27. PROOF OF THEOREM XVIII 
We write (246) in the form 

g{x)=f(x)- I* H{\x-y\)f{y)dy, (249) 

J — oo 

where /(^)==0, a:<0; (/(^) = 0, x>0 (250) 

For x<0,g {x) is understood to be defined by the right-hand side 
of (249). We now set 

{'»')==—= f f(x)€>^dx, y(it)= f (j{x)e^^dx-, 

V27rj-O0 V277 J--0O 

(251) 

in other words, we consider the Fourier transforms 

it), y(s-\- it) 

(regarded as functions of t) of f(x) and g (x) e^'^. 

Since/ (a;) is supposed to satisfy (236), we may state the obvious 

Lemma 1. <f>{u) is holomorphic in the half-plane 5< — a, and 
bounded in every partial half-plane, forinstance, in^ g — (a -}- ^8) /2. 
From (249), x < 0, and from (250), 

|!7(a;)|^J^ \Il{\x-y\)\\f(y)\dy 

<const. \H (\x — y\)\e^ydy. 
For any A between a and 1, we find 

f \H{\x-y\)\e’^<idy<{ \H {\x-y\)\e^y dy 

Jo J -00 

= e^ 

J —00 

y(x) = 0(e-^‘-^i), 


whence, for x<0, 
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A being an arbitrary number less than one. Consequently, the 
Fourier transform has the following property: 

Lemma 2. y {u) is holomorphic in the half-plane 5 > — 1 and 
bounded in every partial half-plane, for instance, in s ^ — (a 4- jS) /2. 

It is the chief point of the following considerations that the 
regularity half-planes of ^ and u have a common ordinate, 
s= — (a4-^)/2. Along such a common ordinate, we have, from 
(249), 

y(M) = ^(M) — ^ e^dx\ H(\x-y\)f{y)dy 

V 27r J —CO J — CO 

= f f{y)dy{ H{\x-y\)e'^dx 

V ZtT J -00 J -00 

= f{y)e«»dy\ H(\t\)e«‘df, 

VZTtJ-oo j —CO 

hence y(u) = (l>(u){l — K (u)}; s= -~(a-l-i8)/2 (252) 


It is permitted to interchange the order of integration along the 
common ordinate, since the double integral converges absolutely 
for 5 = — (a-f-i3)/2. The identity (248) allows us to write (252) in 


the form 


y (M) ^ < f>{u) 
CT^. {U) (U) 


P(u), 


(253) 


5 = — (a-^j8)/2. Here, we know that the left-hand side is holo- 
morphic for 5^ — (a4-)3)/2, while the right-hand side is holo- 
morphic within —(ol + P)I2, Both sides thus define an entire 

function oi u, P (u) being a polynomial of degree 2n, we infer 
from (247) and (247') that this entire function is at most of the 
order of for | u | large. According to a well-known theorem of 
the theory of complex functions it is therefore a polynomial 
Q (u) of degree not greater than n, whence 


(f) (u) = (T_ (u) 


QM 

p(uy 


Q can, according to (247'), not be of degree n since 
being the Fourier transform of a q.s. function, is again a q.s. 
function of t. It is thus proved that (f> {u) has the form (241). 

We must not forget to prove that (241) and (242) actually 
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represent solutions of (228). We now start with (241) and study 
the functions /(a;) defined by (242). (241) and (247') show that 
<f)(u) = 0{\u\~^) for |?^| large, in particular, that ^{8 + it) is a 
q.s. function of t. Taking account of Lemma 1 and of Cauchy’s 
theorem, we see that the abscissa of integration in (242) may 
be moved to the left without change of the integral. For any 
5^ (a-f- j8)/2, is, conversely, the Fourier transform of 

/ (x) Hence and from Parseval’s formula, a = b, A = B, 


j: 


\f(x) \ ^eP‘'*^dx 




+ CO 


I \^dt. 


The right-hand side being bounded for — (a + j8)/2, wo infer, 
on proceeding to the limit s^ — co, that /(:r) vanishes for, at 
least, almost all negative x (in the sense of Lebesgue). Otherwise 
the left-hand side would increase indefinitely as — oo. 

Let us, now, set ^ ^ ^254) 

a-f /S 

2 ’ 

(255) 


2 

and {/ (^) = / — - yM ; s = 

V 2'ITi J H- ioo 


According to (247), y (u) is regular and 0 ( 1 24 1~^) in 5 ^ — (a -f- j3)/2. 
The abscissa of integration can, in (255), be placed arbitrarily 
far to the right, whence we infer that gf (x) vanishes for almost all 
positive X, 

It remains to prove that/(x) and (/(x), found from (242) and 
(255), satisfy (249). In other words, we have to show that, back- 
wards, the relation (252), being satisfied by our and y, implies 
(249). We note that the following functions A, B are the Fourier 
transforms of a, 6, 

a{y) = e‘«'fiy), A {t) = (f>{s + it) = <f>(u), 

b{y)^e^ H{\x-y\), = K'(uj, 

V 27t 


always on the ordinate 5 =— (a + j3)/2. The Parseval relation 
(a,B) = {A, B) gives then, 8= ~(a-f j3)/2, 

^00 2 /•s+ioo 

e*'* H{\x-y\)f{y)dy= <f> {u) k (u) e-^‘^ du. 

Jo v27Tlj8—i<X> 

(256) 


H M P 


6 
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On combining this with (242) and (255), and on taking account 
of (252), we see that (249) and (250), i.e. that (226) is satisfied. 
The right-hand integral in (256) converges absolutely because (f> 
and K are both q.s. functions of t. Both sides in (256) represent, 
therefore, bounded functions of a;. This implies, according to (226), 

The coefficient of the exponent is, here, greater than a. We 
remember, however, that a was subjected to the only restriction 
that no characteristic roots should lie on s= ± a. Our formulae 
would, therefore, not change if a slightly smaller a be employed 
and if jS also bo taken smaller than the original a. There is thus 
no difficulty in establishing (236). 

§28. PROOF OF THE ASYMPTOTIC 
FORM OF THE SOLUTIONS 

According to (241), ^ (i^) can have its poles only among the char- 
acteristic roots. It has, however, at least (n + l) poles within 
1 5 I < a. When, in (242), the abscissa of integration s= — (a -f j8)/2 
is moved to 5 = (a-f j3)/2, the change of the integral amounts to 
— V27T times the sum of the residues of the integrand, contained 
in the strip | 5 | < a. Near a characteristic root u* of order k, the 
integrand is of the form 

a~Jr{u-u*) 

e-'« h (u) {u) ^ («*) 0. 

The residue d^tu^u* has, therefore, the form 

Q* being a polynomial of degree less than k. In the case 
the degree equals precisely A;— 1. Altogether we 

f{x)=f*(x) + r{x), 

where /* (a;) = 0 (e“ '•*’') is a solution of (226'), and where 


1 
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the abscissa of integration now being 5 = (a-f-j3)/2. For a;>0, 
(226') can be written 

/*(a;)-A(/*)^= f H{x+\y\)f*{y)dy, 

J —00 

A being the operator from 0 to +oo, given by (226). This implies, 
according to/^ = 0(e“‘-*^>) and according to (233), 

f* (x)- A {/*),= 0(e~--), 

for cr> 0. The same holds, according to /= A (/), for the function 

r(x)-A(r)^=0(c-“-). (257) 

On the other hand, r {x) 5 = (a -f ^)/2, is q.s., being the Fourier 

transform of along 5 = (a + j3)/2. On applying Schwarz's 

inequality to A (r)^. , the integrand being written 

[U{\x-y\)e-^y]\r{y)e^y], 

we infer, therefore, that r{x)=^0 which completes the 

proof of the corollary. 

§29. NEW PROOF OF THEOREM X 

The essential hypothesis, for the validity of Theorem X, was the 
j)ositivity of the kernel, AT > 0, and (119), which can be written 


/c(0) = 2£f2(0) = l, (258) 

i.e. = 0 is a characteristic root, being at least double. According 
to jff > 0, we have, however, 

k" (0) = 2^^ m {l) dt = (0) > 0 (259) 


which shows that = 0 is precisely double. There are no other real 
roots and, according to | K{it) [ < k (0), no imaginary roots. 
Other complex roots could, however, well exist. The polynomial 
P (u) of (237) can thus be written in the form 

P (u) = u^R (u) S (u), 
where R {u) = — 

8 (u) = {u-{-Ui) 


6-2 
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On making, in (241), the particular choice 

'w^efind 

and, for suitable c, f{x) = x + qoQ-\-o(l) (261) 

(concerning the notation see Theorem VI), since the roots of 
/S (i^) = 0 all lie in the right half-plane 5 > 0, giving rise to residues 
that vanish exponentially as At the same time it is seen 

that no other choice for Q,^_i leads to solutions of the form (261). 


From (239'), (240) and (260), 

U(f)(u)->c; s-^cOy (262) 

and u->0. (263) 

/b(0) 


On the other hand, the first formula (251) is a consequence of 
(242) and holds, here, for all u with 5 < 0, 

(f)(s) = ^_. ( f(x)e^^dx. (264) 

V 277 jo 

In order to determine c, we observe that (264) implies 

— — oc, (265) 

V 27T 

while (261) and (264) imply 

a2<^(s)-> -l=; s->0. (26G) 

V27r 


By combination of (262), (265) and (263), (266), 

''"■‘a:(0)’ 

O’. (0) can be found from (239') and (240), 

(0) . ( - 1 )» exp [ - *] . 


(267) 


(268) 


The abscissa of integration can, here, be moved to the left and 
deformed into another path, consisting of the parts ( — ioo, — ip) 
and {ip, ioo) of the imaginary axis, and of the right half of the 
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circle \ v\=p. The integrals along the first two pieces cancel each 
other, v~^\ogT(v) being an odd function. The expression inside 
of the parenthesis, in (268), equals therefore 

TT 

1 

JogT(pe^^)d^, 

“2 


whence, on proceeding to the limit p = 0, the value 


-ilogT(O) 

follows. From (238), 

T (0) = ( — 


hence 


a_(0) = (-irS(0)lV2H,(0), 


and according to (267), 

It is to be noted that n— 1 is always even, since the conjugate of 
a characteristic root is again such a root. 


Remark on g'oo- Computations which may be left to the reader 
lead to the following value of ^oo hi (261), 


1 (2 . K(it) 


dt. 


In Milne’s case this can be transformed into 


Jo |sin2</> ^ — tan</>j 

The evaluation gives 0*710. 



CHAPTER V 


OTHER PROBLEMS OF 
RADIATIVE EQUILIBRIUM 


§30. PURE ABSORPTION. NON-GRAY MATERIAL 

In this section, a few principal remarks will be made concerning 
the fundamental integral equation in the case of purely absorbing 
material with an absorption coefficient that varies with the wave- 
length. In the case of thermodynamical equilibrium, 

where is a given function of B, 

(269) 

with the following properties, 

^>0; B,(0) = 0; J5,->oo, R->oo (269') 

Only these qualitative properties will be used in the sequel. In 
order to simplify the considerations we suppose that the absorp- 
tion coefficient be of the form 


a,(P) = i3Ha(n (270) 

The material is again supposed to be stratified in parallel 
planes. The variable t defined by 

T = j p(x)oL (a;) dx 

J —00 

is no more the optical depth though being convenient for the 
mathematical discussion. We assume the r-thickness of the slab 
to be finite, r < r^. The true optical depth of a point x, corre- 
sponding to V, is 

T^= p(X.^dx= 

J —00 

writing briefly instead of j8 (v). 
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The equation of transfer is 

cos 0 ^ „ 



while the equation of radiative equilibrium (19) becomes 

= (272) 

The net flux ttF of the radiation of all frequencies is, of course, 
constant. When (271) is multiplied with cos 0 and integrated with 
respect to r and v, we find 

Ji = jFt + const., 

A^=| = ~ f/^cos^edto), (273) 

J 0 '^Pv J 

Avhich represents a generalization of Eddington’s relation (32), 
K having, however, no immediate physical meaning in this case. 

Problem V. Find (r, 6 ) and (t) from (269), (271) and (272) 
when the radiation incident at the surface t = 0 is zero, while the 
radiation incident at the inner face t = t^ has a prescribed value 
independentf of the direction (0 < 7r/2). 

We shall first find the integral equation for B{t). On inte- 
grating (271), 

[//e-^rsec9(’^*->-)-t-(8^sec0 T e-^«’8ec9(t-T) 

IAr,9)=\ ^ 

j^j8„|sec0|J ^274) 

and on integrating, now, through all directions, 

^ = I {pAr-t\)BAB «)} dt + Ov (’•* - r)}. 

(275) 

We introduce the brief notation 

^00 

<!)(*,£) = 1^ p/Ei(P,x)B,{B)dv (27(5) 

and Y{B)=j p,B^{B)dv. (277) 

t This is unessential for the subsequent considerations. 
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When the coefficient of absorption is a known function of depth 
and frequency, O and T are to be considered definite functions of 
their arguments. From (272) and (275), we obtain the funda- 
mental integral equation 

T{B(t)}=J r*<D{|T-«|,J5W}d<+JG(T*-T), 

° (278) 

where G{x)=j (278') 

for the determination of B{t). Once (278) is solved we are able 
to determine = B^{B (r)} and, from (274), (r, 6). 


§31. DISCUSSION OF THE NON-LINEAR 
INTEGRAL EQUATION 


The difficulties connected with the Problem V, pointed out 
already in § 4, compel us to content ourselves with some qualita- 
tive remarks on the solution of (278). Again, it is convenient to 
introduce a brief symbol 


0[B)^^\V ^{\r-tlB{t))dt (279) 

J 0 

for the non-linear integral operator within (278). (278) becomes 

then Y I g = O ( B)^ + (t* - t) (280) 

In the special case of gray material, 1, 0 becomes the linear 
operator L of Schwarzschild, while G {x) simply becomes (x). 

The operator 0 has a fundamental property which may be 
called ‘monotonity’ and which corresponds to positivity in the 
linear case. Since (276) represents, according to (269'), an in- 
creasing function of B, 

Bi(t)^B2(t); Bj^Bj 

0(B^),<0{B,), 

everywhere. We also note that the function T* (B) in (277) is an 
increasing function of B, 


Theorem XIX. Problem V has a solution. The integral equa- 
tion (278) can be solved by successive approximations. 
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Proof. We first study 0{B)^ in the case of a constant B{t). 
For any v, B^ (r) is then constant too. We remember that O (B) 
is obtained from the first term on the right of (275) by multiplica- 
tion with and by integration through the spectrum. This term 
equals, now, o 


whence O (B)^ = Y{B) 




-‘i: 

- ij\ B, E., [^AT*-r)]dv, ...(281) 


in the case of a constant B. We might, together with (276), 
0 =0^, introduce the functions 


(D„(x, ^)=J V-" EJp^x)BAB)dv (276') 

(281) then takes the form 

T (B) = 0 (B),+ 1% (r, B) + JcDg (r^ - T, JS) ( 281 ') 

From (269') we obtain the following propeity of the functions 
O^, in particular of 03 , 

0.2(0;, 5 ) -> 00; B -^00. 

This property shows that, for all large enough constants /?, 

02 ( 0 ;, B)>G(x), 

0<x< or, according to (281'), B= B, that the inequality 

(282) 

0 < T < T^, holds for all sufficiently large constants B. This proves 
the convergence of the successive approximations in (278). 

We set W{B„^,{r)}=0{B,X + W{r*-r) (28.3) 

Bq (t) = 0 . 'F (B) being an increasing function of B, (283) uniquely 
defines a series of functions Bj^{t), 7i — 0 , 1 , 2 , .... We obviously 
have BXT)>B,{r). 

Supposing the inequality B^^ (r) > B^^_^ (r) to be true for one 
particular n, we find from (283), according to the monotonity of 
the operator 0 , 

'F (r)} > O {B„_X + hO (r* - t) = T (5„ (r)), 
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whence (r) > (r). Induction shows thus the general 

validity of this inequality. On the other hand, from Bq (t) < jB, 

'F{5i(T)}<0(5), + 4G(r*-T), 

and according to (282), B^{r)<B. Continuing in this way we 
generally find (t) < B, The limit function 

^(T) = lim B„ (t) 

n=co 

is as before seen to satisfy the integral equation (278), q.e.d. 
The proof of the uniqueness may be omitted. 


§:32. REMARKS ON THE CASE OF 
INFINITE OPTICAL DEPTH 


Theorem I is probably generally true. We shall indicate the proof 
in the present case. 

According to ^ 0, we find as previously 


(^) + i^^sec 6 j B^,dt, (284) 

= radiation from below, 0<7r/2. According to 

= I^(r,r)^e^y^i^(r); 0<7r/2, 


and according to (273), 

K ('r) 


r 

> i^(r)cos^0d(v, 

'^Hv J + 


(285) 


thus yielding 

f I 

^"(t)> j(r,T)cos2 0da^; j= ~-i^{r)dv. 

J Jo Pv 

This is readily seen to imply (r) = 0 for almost all {v, r), because, 
otherwise, K (r) would increase exponentially with increasing r, 
in contradiction to (273). 


Problem VT. Find B and 7^, the incident radiation being 
zero, and the flux constant F ^ 0 being given. 

The integral equation of this problem becomes, according to 
the notation introduced in § 31, 


(286) 
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It is natural to conjecture that this equation has a one*parameter 
family of solutions In the simplest case of gray material, 

j8^=l, (286) becomes Milne’s linear equation. In the present 
general case, the connection between B and F is no longer linear. 
In analogy to the linear case, we could also insert (274), t* = oo, 
directly into the flux integral, whence, according to (276'), 

F = (0} dt (<)} dt (287) 

When this equation is differentiated and account is taken of the 
relations 

- (X, 5) = - <D„ {X, B ) ; Y (B)=% (0, B), 

(286) is again obtained. On inserting (274), t* ~ oo, into (273), we 
find the relation 


J B {()}dt = ^"T-\-const,, 

thasyiolding f 

T--'X> J 0 


This relation might be used in order to determine the parameter 
in B (t, c), when F is given. Differentiation of this relation leads 
again to (287). 


§33. ABSORPTION AND SCATTF:RING. 

SCHWARZSCHILD’S INTEGRAL EQUATION 

When absorption and scattering are simultaneously to be taken 
into account, the determination of the radiation field becomes an 
exceedingly difficult problem. It was pointed out in § 4 that this 
general problem leads to a complicated non-linear integral equa- 
tion of B. In order to derive this equation we must first, by means 
of (16) and (16'), express the Ergiebigkeit in terms of the 
emission = This present section is devoted to this partial 
problem. 

We suppose scattering to be uniform, y= 1. Since the above 
problem concerns only a definite frequency v, we replace, for the 
sake of clearness, the affix i/ by a dash. B^ = B', = .... 
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We introduce the total optical depth 

T= j p{oi a') dx 

J —00 


and set 


A = ~ 


a 4-0- 


1-A = 


a +cr 


In the outermost layers of the sun scattering plays the chief role, 
while in the deeper layers absorption predominates, 

A->1, T->0; A->0, T->oo. (288) 

The incident radiation being zero, we obtain, on integrating 
(16) and on inserting the intensity into (16'), = 

Schwarzschild’s integral equation 

J' (r) = A (r) A {JX + {1 - A (r)} B' (r) (289) 

for the determination of the Ergiebigkeit, A being Milne’s 
operator (52). Milne has shown that an equation of the type of 
(289), together with (288), holds under much more general con- 
ditions than under local thermodynamical equilibrium. We con- 
tent ourselves with a few qualitative remarks on the solution. 
In order that (289) have a i)ositive solution, the A-solution must 
be finite. We prove that, in general, the Neumann series con- 
verges. 


Theorem XX. Under the hypothesis (288), the N-seriea of 
(289) converges if ^ q 

holds for some a < 1 . 


Proof, First we study the homogeneous equation corresponding 
to (289), but with a constant A = A. Since (230) represents the 
characteristic function of Milne’s equation /= A(/), 



1 -i-?4 

l^u 


will be the characteristic function of 

/(t)=AA(/),. (290) 

We know that, for A= 1, = 0 is the only characteristic root in 

the strip 1 5 | < 1, and that, for A < 1, there are precisely two roots 
in that strip, being real and of opposite sign. When u increases 
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from 0 to 1, increases from A to +oo. For an arbitrarily 
given positive j8< 1, we can, therefore, find a A< 1 such that p 
becomes a characteristic root. According to the theory of Chapter 
IV, there exists a solution of (290), satisfying 

f{r)-ePr, (291) 

for T large. The following choice of /8 is convenient for our pur- 

a<jS<l, (291') 

a referring to the hypothesis of the theorem. 

From (290), 

/+ c — AA (/+ c) = ^ 1 — ^ j/+ c (1 — AA ( 1 )}, . . . (292) 

c being a constant. We observe that A = A (r) is always less than 
one and tends to zero as t->oo. According to the hypothesis 
made about -B'(t), and according to (291) and (291'), the first 
term on the right of (292) will, for all large r, be greater than 
{1 — A(t)}jB'(t). Onthcotherhand we can, according to a (1)< 1, 
choose c so large that, in the remaining finite r-range, the right- 
hand side of (292) becomes greater than (1 — A) B\ This choice of 
c leads thus, for ail r, to 

J (t) > a (r) a ( j), + {1 - A (r)} jB' (r), (293) 

J J is, of course, positive when c is sufficiently large. 

This proves, in the usual way, the convergence of the A-series. 
Denoting the partial sums of that series by (r), Jq = 0, we have 

e/,^i = AA(JJ-f(l~A)i?'. 

On subtracting this from (293), 

The positivity of A implies thus Jn<J for all n and for all r. 
From Jq<Ji<J 2 <^->, the existence of the limit function «/', 
solving (289), is inferred, q.e.d. 

It may be mentioned that the homogeneous equation resulting 
from (289) has no solution (except zero) provided that A(t)-> 0 
converges sufficiently rapidly as t->oo. 
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§34. MILNE’S MODEL OF A PLANETARY NEBULA 
IN RADIATIVE EQUILIBRIUM 

There is another problem of radiative equilibrium worth men- 
tioning because of different boundary conditions. A spherical 
gaseous shell is illuminated by a source of light located in the 
centre of the shell. What is the distribution of light in the shell 
when in radiative equilibrium? We confine our attention to the 
case of purely absorbing gray material, or to the formally equi- 
valent case of monochromatic radiative equilibrium. Let, for 
the sake of simplicity, the shell be infinitesimally thin, the 
optical thickness being kept finite (see the beginning of §5), 
i.e. we neglect curvature. 

We introduce the same variables r and 6 as in the previous 
models of a stellar atmosphere. The boundary conditions are 
then as follows : 

(a) There is no radiation incident on the outer face, t = 0, 

/(O,0) = O; e>7T/2, 

(b) The inner face receives the normally incident parallel 

radiation of the point source 0, the net flux, at t = being ttS, 

(c) This, however, does not ex- 
haust the radiation incident at the 
inner face. It also receives radia- 
tion emergent from other points 
of the inner face. The radiation 
not being weakened while travel- 
ling through empty space, the 
figure shows that w^e have 

I (t*, 6) = I {T'^,7T — d); 0<7 tI2 
at the inner face. 

This apparently represents a new tyi)e of boundary condition. 
It might be regarded as a special case of 

I (r^, 0) = / (r*, 77 - 0) + i (0) ; 0 < 7r/2, (294) 

i (0), the excess of the incident over the emergent radiation, being 
given. When the finite size of the central star in O is taken into 
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account, i (6) represents a finite function. The case of parallel 
radiation normally incident at t = t* must, as in § 17, be treated 
as a limit case.! 

We now derive Milne’s integral equation of the problem, 
under the boundary conditions (a) and (294). Schwarzschild’s 
fundamental integral equation of the plane model, with the 
boundary condition (a), and with any radiation incident at r = r*, 
was given by (64) and (64'). According to (294), the last term in 
(64') splits into the summands 


4itJ+ 477 J + 

(295) 

The radiation in the first term comes through the shell from its 
outer parts, 

/(T^7r-0) = sec0 r 0<7r/2. 

Jo 

Inserting this into (295) and using (49), n=l, we obtain the 
integral equation 

= (296) 

477 J 4 


where 


J(t)[E(\r^t\) + E(2T*-T-t)]dL 
Jo 


.(296') 


On proceeding, as in the beginning of § 17, to the limit case of 
normally incident radiation i of net flux tt/S, we obtain Milne’s 

equation o 

J(r) = L*(J), + ~er-r* (297) 


The net flux ttF, being constant in the shell, is found from (294), 

F = - f i(i/)cos0dco. 

^J-h 

In the present limit case we find F = S, 


t For a finite central star, (294) does not precisely represent the condition of 
the problem, because the central star covers part of the shell. In the limit case, 
the condition is, of course, exact. It will, however, be of advantage in the sequel 
to consider the mathematically more general condition (294). 
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It might be mentioned that the integral equation can easily 
be brought into the simpler form 

l‘2r* S 

^(r) = iJo J(t)E(\r-t\)dt + -e-^^*-^K 

J (t) being defined in the larger interval (0, 2t*) by reflection at 
T — T*, This shows, as before, that the solution is unique and given 
by the convergent iV-series. The same is, of course, true in the 
more general case (296), since this equation admits of an analogous 
transformation . 

The radiation emergent from the outer face is given by the 
usual formula (47), r = 0. It is of interest to know how the model, 
S (more generally i{d)) being given once for all, behaves as the 
shell becomes more and more opaque (t*->oo). Milne's approxi- 
mate formulae suggest that the distribution of light, in particular 
of the emergent radiation, becomes the same as in model Ja, 
F = 8. The mere fact that the flux F= 8, in the shell, remains the 
same for all r*, shows already that the shell cannot become dark 
as r^->oo. The central star is, of course, seen through the shell, 
but its light is dimmed by absorption (or scattering). The greater 
part of the light of the shell comes, therefore, from the shell 
itself. We now give a rigorous proof of the above-mentioned fact. 

Theorem XXI. When t*->oo, the solution of (297) ap- 
proaches the function -8f{T) 

i.e. the Ergicbigkeit in the case la, F = 8, the convergence being 
uniform in every finite r-interval. 

Proof, It is convenient to split the equation (297) into two 


simultaneous integral equations, 

j7(T) = L(?),+|e--*. (298) 

L being Schwarzschild’s operator (64), and 

h ( t) = L* {h)^ + -^ 9 (0 E (2t* -T-t)dt, (298') 

L* being the present operator (296'). We obviously have 

= ^ + (299) 
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I 7 , h depend, of course, upon r* too. The physical meaning of (298) 
is that (] represents the Ergiebigkeit when the radiation received 
from the other portions of the inner face is neglected. Finally, 
(298') is the equation of the iHoblem, where the true boundary 
condition (294) is taken account of, i (6) being identified with the 
radiation emergent in the former case. 

According to the symmetry property (210) of Schwarzschild’s 
operator, (j(t) = (J (r* — r) becomes the solution of 

{/ (r) = L ((/), + (300) 

If in the A-series representing {/ all integrals are extended up to 
infinity instead of to all terms obviously increase, i.c. (j (r) 
is smaller than the A-solution of the same equation, but with 
Milne’s operator A. The latter solution has, however, in § 17 been 
shown to be bounded, {/(r) is thus uniformly bounded for all r*. 
Therefore, from (298'), 

A (t) < L* {hi + 1 (r* - r), (301) 

c> 0 being a suitable constant. 

We now make a crude comparison with model la. The theory 
of model la shows that there 

I (t*, 0) — I (r*, 77 — 0) > c' cos 0, 

6 < 77 / 2 , c' being a suitable i)ositive constant. Inserting this into 
(296), J being identified with the function /(r) of Problem la, 
we get 

/(t)>L*(/),+ -io’3(T*-T) (302) 

(301) implies that }i(j) lies below the A-series corresponding to 
the last term in (301). Similarly, /(r) lies above the A-serics 
corresponding to the last term in (301). From Ag < 2 A 3 we find, 
therefore, ^ 

^ ('^) < ; 7 / 

C/ 

thus yielding, altogether, an inequality 

J(T)<AS^a(T+l), (303) 

a being independent of t*. 


H MP 


7 
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(303) will now be used to show the uniform continuity of J (t) 
for all r*. Let us first fix an (otherwise arbitrary) interval 
We choose T*>k and express, according to (290) and 
(297), the difference 

J(t")-J(t'), 

in terms of integrals. Arranging these integrals in a convenient 
way, we find 

I ./ (t") - ./ (r') 1 < " J >, + 1 ) I A’ ( 1 t" - I ) - ( I r' - i I ) 1 dt 

+ o r (<+l)| A(2T*-T''-0-7i!;(2T*-r'-<)|d< 

The second term is readily seen to be smaller than 

The above inccpiality signifies, therefore, that in an arbitrarily 
given T-interval (0, ^), the functions J (r) are uniformly con- 
tinuous, for all T*>k+l, 

According to Arzela^s selection theorem, any sequence of 
functions J (t,t*), taken for an arbitrary sequence of numbers 
r^*->co, possesses a subsccpicnce that converges uniformly in 
every finite r-interval. Let J (r) be the limit function of such a 
subsequence. It is on account of the uniform inequality (303), 
that, in (297), the integration and limit process can be inter- 
changed, thus yielding 

J(r) = A(j\, J^O, 

A being Milne’s operator. Theorem V show.s, therefore, that J ( t ) 
is a constant multiple of /(t). The constant of proportionality is 
easily determined by means of the fact that the flux constant is 
always jS, On applying the auxiliary theorem of § 13 to (297) and 
to the relation 

l = L^(ll-h iJS, (r) + (2r* - r), 

we find 2 J (t) i^) + (2t* - t)] = /S ( 1 - e-^*). 
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According to (303), we may again proceed to the limit t*->oo 
under the integral sign, 

2jV(T)i?2(T) = S, 

thus yielding (T) = lSf(r), Since, thus, all converging sub- 
sequences of functions J, t^->oo, have the same limit function, 
we infer that J (r, r^) converges towards the mentioned limit 
function as the convergence being uniform in every finite 

T-interval, q.e.d. 

When the radiation from other portions of the inner face is 
neglected, i.e. when (6) is the only boundary condition at t = t^, 
the integral equation is (298), ^(r) being the Ergiebigkeit. It is, 
in this case, physically obvious that the shell becomes com- 
pletely dark as t*->oo, i.e. that 

f/^0, r*->oo. (304) 

This can be rigorously proved as follows. Applying the auxiliary 
theorem of § 13 to (208) and (300), we obtain 

Knowing from §22 that m< 1 and that as t*->oo, we 

infer that the last integral in this equation tends to zero as 

T*->00. 


In (298), the last term tends to zero as As to the 

integral representing L (g), we split the integrand into the two 
factors / E(\T-t\)\ 




VE^(t) 

By application of Schwarz’s inequality, 


L {gV < i jj it) dt. 

The first factor has been shown to tend to zero as t*->oo. Ac- 
cording to the boundedness of g, the second factor is less than a 
constant times 


(298) shows thus that (304) is true. 


i: 


E^it) 


dt<oo. 


7-2 
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It should be emphasized that this bibliography refers only to 
mathematical contributions to the problems in question, and to 
the original memoirs concerning the astronomical and geophysical 
applications. A comprehensive exposition of the theory of radia- 
tive equilibrium and of the approximate solution of the problems 
is given in Milne (1). 

To §§ 1-4. The theoretical foundations of the theory of radi- 
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ative equilibrium are due to Schwarzschild (1), (2). He made the 
first fundamental applications to the sun’s atmosphere (law of 
darkening on the sun’s disk, origin of the Fraunhofer lines in the 
spectrum). To the interior of a star, the theory was first applied 
by Eddington. The mathematical problems arising from the 
stellar interior are, however, not taken up in this tract. 

It is necessary to say a few words about the customary intro- 
duction of the fundamental quantities of the radiation field. The 
use of several differentials ds, da^ doj, dm, dr, especially in the 
definition of intensity, is rather troublesome. From the point of 
view of the mathematician, it lacks both rigour and beauty. The 
entire subject ought to be treated anew, and the fundamental 
equations be derived in a rigorous way, the main tool being the 
general theory of measure (in the sense of Lebesgue-Radon). 

To § 5. The purely absorbing gray model was first introduced 
by Schwarzschild (1), while the model of a purely scattering 
atmosphere is due to Schuster. The radiation field was originally 
replaced by two antiparallel streams (Schuster-Schwarzschild- 
approximation). The Schuster model was rigorously set up and 
treated in Schwarzschild’s fundamental memoir (2) (case of 
uniform scattering). Finally, the rigorous equations for the ab- 
sorbing gray atmosphere were approximately solved by Milne, 
who also studied the spectral distribution of the emergent light 
(cf. Milne (1)). 

To § 7. The reduction to integral equations of boundary value 
problems of the elementary theory of radiation is due to Hilbert, 
who used it in his attempt to prove Kirchhoff ’s laws. The funda- 
mental integral equation (289) for the Ergiebigkeit, obtained from 
(16) and (16') in the plane case, is due to Schwarzschild (2). 
The integral equations (53) and (54) of Problems 1, II have been 
given by Milne (1), (2). Milne established (53) originally in a 
different form, which is obtained easily by partial integration 
of the integral A (J)^. On writing 

C(r)=j^J{t)dt, 


(305) 
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we find 

^ (306) 

The importance of the positivity of the kernels has already 
been realized by Hilbert and Schwarzschild. 

To § 9. This method was first given in the author’s paper (1). 
The existence of a solution could, of course, also be obtained from 
Schwarzschild’s theorem about the solution of (64'), the last 
term being (66). Schwarzschild proved that the solution tends 
towards a limit function as t*->oo. Since (cf. Theorem XVI) 
for every the solution lies below r + 1 and above r, it could 
easily be seen that limit process, t^->oo, and integration can be 
interchanged in Schwarzschild’s integral equation. The limit 
function lies thus between r and r+l and satisfies Milne’s 
equation (cf. Kostitzin, Hopf (1)). In view of the central 
significance of Model la, however, a direct treatment seems 
preferable. The relation c — \F was independently proved by 
Bronstein (1) and by the author (6). The simpler proof given here 
is the author’s. 

To § 10. The uniqueness was proved by the author (2). 

To §§11, 12. Concerning the more general equation, treated 
there, see also Hopf (3). 

To §14. Formula (154) for the boundary temperature was 
independently proved by Bronstein (2) and by the author (6). 
The common root of both proofs is the general formula (151). 
This formula is equivalent with the following fact, proved by 
Ambarzumian in the case of Milne’s equation. The resolvent 
kernel K (t, t) of the given kernel (| t — ^ |) satisfies the rela- 
tion K {0,t)=f (t)/fQ. Concerning Theorem X see the author’s 
paper (8). 

To §§ 15, 16. The increasing of q (r) was proved by the author 
[7), The equation (176) is due to Bronstein (3). The numerical 
-esults show that, at least for small and large t, Eddington’s 
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approximation of ^(t) is the most accurate one (cf. Milne (1), 
Table I). 

To §17. Model 16, for parallel incident radiation, was con- 
sidered by Milne and Eddington. The case == 0 applies, according 
to Milne (2), to the upper air in the earth’s atmosphere, while the 
case -F > 0 is realized in the case of close binary stars (reflection 
effect). For the approximate solution see Milne (1). 

To §21. Concerning Arzela’s selection theorem, used also in 
§ (35), see Courant-Hilbert, Methoden der Matlmnatischen Physik, 
Berlin. 

To §§21-23. Theorem XVI is, in the case of uniform scat- 
tering, due to Schwarzschild (2), as well as the convergence 
towards a limit model, expressed in Theorem XVII. The 
Schuster-Schwarzschild approximate formula is 


i.e. the mean of Schwarzschild’s limits. Schwarzschild computed 
the correction on replacing the integral equation by a system of 
linear algebraic equations. The approximate expression for J (r), 
given in (217'), can numerically be computed by using a very 
good approximation of ^(r), for instance Eddington’s (Milne (1)), 
which probably differs less than 0*5 per cent, from the true q (r). 
A comprehensive account of the Schuster-Schwarzschild pro- 
blem is given in Milne (1), § 14. 

To §§ 24-30. The solution by means of Fourier-Laplace in- 
tegrals was given by Norbert Wiener and the author. The 
explicit, though very complex, formula for the law of darkening, 
resulting from this theory, is of main interest in astrophysics. It 
must, however, be mentioned that this method, though solving 
Problem la explicitly, has but a limited applicability, as concerns 
many other problems of radiative equilibrium. The comparison 
methods, resting on the positivity of the kernel, furnish a simpler 
and more general access to a physically satisfactory (though not 
explicit) solution. 



104 BIBLIOGRAPHY AND SUPPLEMENTARY REMARKS 
The simpler equation (226'), 

B(t)E{\T^t\)dt, (307) 

J —00 

can be interpreted as the approximate form at great depth of 
Milne’s exact equation. The only characteristic roots being u = (^ 
(double), (r) = ar + 6 is suggested to be the only solutions of the 
type (236). On using (305), J = B, the above equation becomes 

c (t) = P ^ e-idt (307') 

Jo 

This equation has been studied by Littlewood and by Hardy 
and Titchmarsh. An account of their results is given in Milne 
(1). The latter two authors have proved that the only continu- 
ously differentiable solutions of type (236) are given by 

C (r) = ar^ -i-br + C, 

The more general equation, obtained from (307') on replacing c" 
by an arbitrary function vanishing like an exponential function 
as <-> 00 , has been treated by the author (4). 

The proof of Theorem X, § 29, was given by the author (8). 

To § 34. About the relation between radiation and tempera- 
ture, in absence of local thermodynamical equilibrium, see Milne 
(1), § 19! Our A is called there 1/1 +r]. 

To § 35. The ‘ planetary nebula’ problem was treated in Milne 
(3). Reference to the former work by Jeans and Gerasimovich 
is found at the end of this paper. 



LIST OF SOME FORMULAE 
OF PHYSICAL INTEREST 


Pure absorption. Gray material in strict radiative and local 
thermodynamical equilibrium. Model la: 

^ ^ = P’ {t + ? (t)}, ^ g (t) < (Jao, 

7tF being the net flux, being the temperature at the optical 
depth T. 

T ^ T ^ 

0 ~ 4 ^ <5 j 

Tq being the surface, the effective temperature. 

_ 1 ^ \ . 

’’■Jo — 

The function 

V» (») = (1 + 1) exp j f 'liSHl - j 

can once for all be tabulated. The law of darkening is (see end of 
§25) /(U,0) = |F<I>(sec^). 

This holds also in the absence of local thermodynamical equili- 
brium. 

Model I b : Parallel incident radiation of normal net flux ttS, 
the angle of incidence being 6'. The absorption coefficient for the 
incident radiation is supposed to be a constant fraction n of the 
general absorption coefficient. The following formulae give the 

upper and inner limiting temperature Tq , 2’a, ^ R = ^ T*^ , 

Bq = ^ nS^ (n seed'); Ba^ = ^S cos d'O {n sec 6') 

BJB^ = n sec d'iVs. 


and 
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